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Abstract 

We calculate {g — 2) of the muon and the QED coupling a(M^), by improving the 
determination of the hadronic vacuum polarization contributions and their uncertainties. 
We include the recently re-analysed CMD-2 data on e^e~ vt+tt^. We carefully combine 
a wide variety of data for the e^e~ production of hadrons, and obtain the optimum form 
of R{s) = <7[|jj^(s)/(Tpt(s), together with its uncertainty. Our results for the hadronic 
contributions to 5 - 2 of the muon are a)^^^'^ = (692.4 ± 5.9exp ± 2.4i.ad) x 10"^° and 
^had.NLO ^ (_9 8 ^ o.lexp ± O.Orad) X IQ-^^, and for the QED coupling Aa[f2d(^l) = 
(275.5±1.9exp±1.3rad) X 10""^. These yield {g-2)/2 = 0.00116591763(74), which is about 
2.4cj below the present world average measurement, and a(M|)-i = 128.954 ± 0.031. We 
compare our [g — 2) value with other predictions and, in particular, make a detailed 
comparison with the latest determination of {g — 2) by Davier et al. 
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1 Introduction 

Hadronic vacuum polarization effects play a key role in the prediction of many physical quan- 
tities. Here we are concerned with their effect on the prediction of the anomalous magnetic 
moment of the muon, = {g^ — 2)/2, and on the running of the QED coupling to the Z boson 
mass. We explain below why it is crucial to predict these two quantities as precisely as possible 
in order to test the Standard Model and to probe New Physics. 

First, we recall that the anomalous magnetic moments of the electron and muon are two 
of the most accurately measured quantities in particle physics. Indeed the anomalous moment 
of the electron has been measured to a few parts per billion and is found to be completely 
described by quantum electrodynamics. This is the most precisely tested agreement between 
experiment and quantum field theory. On the other hand, since the muon is some 200 times 
heavier than the electron, its moment is sensitive to small-distance strong and weak interaction 
effects, and therefore depends on all aspects of the Standard Model. The world average of the 
existing measurements of the anomalous magnetic moment of the muon is 

a^'P = 11659203(8) x 10"^°, (1) 

which is dominated by the recent value obtained by the Muon g — 2 collaboration at Brookhaven 
National Laboratory pp. Again, the extremely accurate measurement offers a stringent test of 
theory, but this time of the whole Standard Model. If a statistically significant deviation, 
no matter how tiny, can be definitively established between the measured value a°^P and the 
Standard Model prediction, then it will herald the existence of new physics beyond the Standard 
Model. In particular the comparison offers valuable constraints on possible contributions from 
SUSY particles. 

The other quantity, the QED coupling at the Z boson mass, Mz, is equally important. 
It is the least well known of the three parameters (the Fermi constant G^, Mz and a{Mz)), 
which are usually taken to define the electroweak part of the Standard Model. Its uncertainty 
is therefore one of the major limiting factors for precision electroweak physics. It limits, for 
example, the accuracy of the indirect estimate of the Higgs mass in the Standard Model. 
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The hadronic contributions to (? — 2 of the muon and to the running of a{s) can be calculated 
from perturbative QCD (pQCD) only for energies well above the heavy flavour thresholds^. To 
calculate the important non-perturbative contributions from the low energy hadronic vacuum 
polarization insertions in the photon propagator we use the measured total cross section^ 

Kadis) = atot(e^e" ^ 7* ^ hadrons), (2) 

where the superscript is to indicate that we take the bare cross section with no initial state 
radiative or vacuum polarization corrections, but with final state radiative corrections. Alter- 
natively we may use 

R{s) = ^ , (3) 

where a^t = 4'7ra^/3s with a = a{0). Analyticity and the optical theorem then yield the 
dispersion relations 

ar'^^ = (^y Tds^^^, (4) 





roc 


3tc J J 






roo 






3n J 





for the hadronic contributions to = {g^ — 2)/2 and Aa{s) = 1 — a/a{s), respectively. The 
superscript LO on denotes the leading-order hadronic contribution. There are also sizeable 
next-to-leading order (NLO) vacuum polarization and so-called "light-by-light" hadronic con- 
tributions to a^, which we will introduce later. The kernel K{s) in is a known function (see 
(|45|) ) . which increases monotonically from 0.40 at m^o (the ir^j threshold) to 0.63 at s = 4m^ 
(the 7r+7r~ threshold), and then to 1 as s — > 00. As compared to (0) evaluated at s = M|, 
we see that the integral in is much more dominated by contributions from the low energy 
domain. 

At present, the accuracy to which these hadronic corrections can be calculated is the limiting 
factor in the precision to which g — 2 of the muon and a(M|) can be calculated. The hadronic 
corrections in turn rely on the accuracy to which R{s) can be determined from the experimental 
data, particularly in the low energy domain. For a precision analysis, the reliance on the 
experimental values of R{s) or o-^s,d{s) poses several problems: 

• First, we must study how the data have been corrected for radiative effects. For example, 
to express R{s) in (@)) and (0) in terms of the observed hadron production cross section, 
crhad(s), we have 

^ ' AnaySs \a{s) J AnaySs ' ^ ^ 

if the data have not been corrected for vacuum polarization effects. The radiative correc- 
tion factors, such as (a/a(s))^ in (jH)), depend on each experiment, and we discuss them 
in detail in Section |21 



^In some previous analyses pQCD has been used in certain regions between the flavour thresholds. With 
the recent data, we find that the pQCD and data driven numbers are in agreement and not much more can be 
gained by using pQCD in a wider range. 

^Strictly speaking we are dealing with a fully inclusive cross section which includes final state radiation, 
e+e" hadrons(+7). 
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• Second, below about ^/s ~ 1.5 GeV, inclusive measurements of cr°g^d(s) are not available, 
and instead a sum of the measurements of exclusive processes {e~^e~ tt+tt^, tt^tt'tc^, 
K^K~ , . . .) is used. 

• To obtain the most reliable 'experimental' values for R{s) or 0"{Jad('5) "we have to combine 
carefully, in a consistent way, data from a variety of experiments of differing precision and 
covering different energy intervals. In Section El we show how this is accomplished using 
a clustering method which minimizes a non-linear function. 

• In the region 1.5 < ^/s < 2 GeV where both inclusive and exclusive experimental deter- 
minations of (Jhadi^) have been made, there appears to be some difference in the values. 
In Section El we introduce QCD sum rules explicitly designed to resolve this discrepancy. 

• Finally, we have to decide whether to use the indirect information on e^e^ hadrons ob- 
tained for ^/s < rrir, via the Conserved- Vector- Current (CVC) hypothesis, from precision 
data for the hadronic decays of r leptons. However, recent experiments at Novosibirsk 
have significantly improved the accuracy of the measurements of the e+e~ —>■ hadronic 
channels, and reveal a sizeable discrepancy with the CVC prediction from the r data; 
see the careful study of [2 . Even with the re-analysed CMD-2 data the discrepancy 
still remains [3]. This suggests that the understanding of the CVC hypothesis may be 
inadequate at the desired level of precision. It is also possible that the discrepancy is com- 
ing from the e~^e~ or r spectral function data itself, e.g. from some not yet understood 
systematic effect.^ 

The experimental discrepancy may be clarified by measurements of the radiative return'* 
events, that is e~^e~ — > 7r+7r~7, at DA$NE [B] and BaBar jH]. Indeed the preliminary 
measurements of the pion form factor by the KLOE Collaboration P compare well with 
the recent precise CMD-2 tt+tt^ data pHl ITT] in the energy region above 0.7 GeV, and 
are significantly below the values obtained, via CVC, from r decays [2]. We therefore do 
not include the r data in our analysis. 

We have previously published ^ short summary of our evaluation of (^, which gave 

^had,LO ^ ^g83_^ ^ 5_g^^^ ^ 2.0rad) X 10"^°. (7) 

When this was combined with the other contributions to (7 — 2 we found that 

a™ = {g- 2)/2 = (11659166.9 ± 7.4) x 10"^°, (8) 

in the Standard Model, which is about three standard deviations below the measured value 
given in (jT)). The purpose of this paper is threefold. First, to describe our method of analysis in 
detail, and to make a careful comparison with the contemporary evaluation of Ref. Second, 

■^The energy dependence of the discrepancy between e^e~ and t data is displayed in Fig. 2 of 3 . One 
possible origin would be an unexpectedly large mass difference between charged and neutral p mesons, see, for 
example, 0]. 

*See for a theoretical discussion of the application of "radiative return" to measure the cross sections for 
e+e~ TTTT, KK, ... at and B factories ^\7\. 
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the recent CMD-2 data for the e+e" ^ Tr+Tr", -k+ti-ti^ and K^Kl channels [HIIISIIII has just 
been re-analysed, and the measured values re-adjusted ^U]. We therefore recompute a^^'^'^'^ to 
see how the values given in ((7j) and (jSj) are changed. Third, we use our knowledge of the data 
for R[s) to give an updated determination of Aahad(s), and hence of the QED coupling a;(M|). 

The outhne of the paper is as follows. As mentioned above. Section |2l describes how to 
process and combine the data, from a wide variety of different experiments, so as to give the 
optimum form of -R(s), defined in Q. In Section |21 we describe how we evaluate dispersion 
relations (jlj) and (0), for a}^'^'^'^ and Aohad respectively, and, in particular, give Tables and plots 
to show which energy intervals give the dominant contributions and dominant uncertainties. 
Section m shows how QCD sum rules may be used to resolve discrepancies between the inclusive 
and exclusive measurements of R{s). Sectional contains a comparison with other predictions of 
(7—2, and in particular a contribution- by-contribution comparison with the very recent DEHZ 03 
determination In Section IHl we calculate the internaf hadronic light-by- light contributions 
to a^. Section [7| describes an updated calculation of the NLO hadronic contribution, a}'^'^'^^'^ . 
In this Section we give our prediction for (7 — 2 of the muon. Section |H1 is devoted to the 
computation of the value of the QED coupling at the Z boson mass, a(M|); comparison is 
made with earlier determinations. We also give the implications of the updated value for 
the estimate of the Standard Model Higgs boson mass. Finally in Section IHl we present our 
conclusions. 



2 Processing the data for e+e hadrons 

The data that are used in this analysis for -R(s), in order to evaluate dispersion relations 
^ and (0), are summarized in Tabled for both the individual exclusive channels (e^e~ — >■ 
TT+TT", Tr'^TT^vr^, K'^ K~ , . . .) and the inclusive process {e~^e~ — 7* — hadrons)^. In Sections 12.11 - 
12.31 we discuss the radiative corrections to the individual data sets, and then in Section 1231 we 
address the problem of combining different data sets for a given channel. 

Incidentally, we need to assume that initial state radiative corrections (which are described 
by pure QED) have been properly accounted for in all experiments. We note that the interfer- 
ence between initial and final state radiation cancels out in the total cross section. 

2.1 Vacuum polarization corrections 

The observed cross sections in e~^e~ annihilation contain effects from the s-channel photon 
vacuum polarization (VP) corrections. Their net effect can be expressed by replacing the QED 
coupling constant by the running effective coupling as follows: 

a{sY. (9) 

On the other hand, the hadronic cross section which enters the dispersion integral representa- 
tions of the vacuum polarization contribution in (jH) and (0) should be the bare cross section. 

^In this notation, the famihar hght-by-hght contributions are called external; see Section |S1 
complete compilation of these data can be found in 
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We therefore need to multiply the experimental data by the factor 

if no VP corrections have been applied to the data and if the luminosity is measured correctly by 
taking into account all the VP corrections to the processes used for the luminosity measurement. 
These two conditions are met only for some recent data. 

In some early experiments (DM2, NA7), the muon-pair production process is used as the 
normalization cross section, cxnorm- For these measurements, all the corrections to the photon 
propagator cancel out exactly, and the correction factor is unity: 

ap = cfp = 1- (11) 

However, most experiments use Bhabha scattering as the normalization (or luminosity-defining) 
process. If no VP correction has been applied to this normalization cross section, the correction 
is dominated by the contribution to the t channel photon exchange amplitudes at tram, since 
the Bhabha scattering cross section behaves as da/dt oc a^/t^ at small |t|. Thus we may 
approximate the correction factor for the Bhabha scattering cross section by 

^ («(tmin))'. (12) 

In this case, the cross section should be multiplied by the factor 

(a/a(tmin))^ 




ap = = ,r;:u'\.. = I ^ 1 ' (13) 



where 

tmin = -S(l - COS^cut)/2. (14) 

If, for example, | cos^cutl — 1? then a(tmin) — and the correction factor (fT^ would be nearer 
to (fTUj) . On the other hand, if | cos6'cut| ^ 0.5, then Q!(tmin) ~ ^(■s), and the correction (fT^ 
would be near to (|TT| . 

In most of the old data, the leptonic (electron and muon) contribution to the photon vacuum 
polarization function has been accounted for in the analysis. [This does not affect data which 
use a{fi^fi~) as the normalization cross section, since the correction cancels out, and so (fTT|l 
still applies.] However, for those experiments which use Bhabha scattering to normalize the 
data, the correction factor (fT^ should be modified to 



(«Ktmin)/«(t„n„))^ ' ^ ' 



where ai{s) is the running QED coupling with only the electron and muon contributions to the 
photon vacuum polarization function included. In the case of the older inclusive R data, only 
the electron contribution has been taken into account, and we take only / = e in ()15|): 



{ae{s)/a{s)f 
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We summarize the information we use for the vacuum polarization corrections in Table |21 
where we partly use information given in Table III of [221 and in addition give corrections for 
further data sets and recent experiments not covered there. It is important to note that the 
most recent data from CMD-2 for 7i~^tt~, tt+tt^tt'^ and KgR^^, as re-analysed in ^U], and the 
KgKl data above the [IH], are already presented as undressed cross sections, and hence are 
not further corrected by us. The same applies to the inclusive R measurements from BES, 
CLEO, LENA and Crystal Ball. In the last column of Table |21 we present the ranges of vacuum 
polarization correction factors Cvp, if we approximate - as done in many analyses - the required 
time-like a{s) by the smooth space-like a{—s). The numbers result from applying formulae ()1U|). 
(fTT|) . (fT3j). (Uni) as specified in the second to last column, over the energy ranges relevant for the 
respective data sets''. The correction factors obtained in this way are very close to, but below, 

°n. 1.1 
U 

1.075 



1.05 



1.025 



1 



0.975 



0.95 



0.925 
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Figure 1: Vacuum polarization correction factor in the low energy regime. The continuous 
line is the full result as applied in our analysis, whereas the dashed line is obtained when using 
a{—s) as an approximation for a{s). Both curves are for cos6'cut = 0.5 whereas the dotted hues 
are obtained for cos^^ut =0.8. 

^To obtain these numbers we have used the parametrization of Burkhardt and Pietrzyk for a(g^) in 
the space- hke region, < 0. 
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one, decrease with increasing energy, and are very similar to the corrections factors as given in 
Table III of jHS]- However, for our actual analysis we make use of a recent parametrization of 
a, which is also available in the time-like regime |l(Jlj . For the low energies around the uo and 
resonances relevant here, the running of a exhibits a striking energy dependence, and so do 
our correction factors Cvp. We therefore do not include them in Table El but display the energy 
dependent factor in Fig. [T] For comparison, the correction using space-like a, a(— s), is 
displayed as dashed and dotted lines for cos6'cut = 0.5 and 0.8 respectively. 

For all exclusive data sets not mentioned in Table El no corrections are applied. In most 
of these cases the possible effect is very small compared to the large systematic errors or even 
included already in the error estimates of the experiments. For all inclusive data sets not cited 
in Table |21 (but used in our analysis as indicated in Table we assume, in line with earlier 
analyses, that only electronic VP corrections have been applied to the quoted hadronic cross 
section values. We therefore do correct for missing leptonic (yU, r) and hadronic contributions, 
using a variant of (fTUI) without the electronic corrections: 

{a/ae{s)Y 




ap = = = 1 5^ 1 ■ (17) 



This may, as is clear from the discussion above, lead to an overcorrection due to a possible 
cancellation between corrections to the luminosity defining and hadronic cross sections, in 
which case either (if (Xnorm = cfnn) or (if cinorm = o"ee) should be used. However, those 
corrections turn out to be small compared to the error in the corresponding energy regimes. In 
addition, we conservatively include these uncertainties in the estimate of an extra error 6a^^, 
as discussed below. 

The application of the strongly energy dependent VP corrections leads to shifts Aa^^ of 
the contributions to as displayed in Table El Note that these VP corrections are significant 
and of the order of the experimental error in these channels. In view of this, the large positive 
shift for the leading tt+tt" channel — expected from the correction factor as displayed in Fig.^ 
— is still comparably small. This is due to the dominant role of the CMD-2 data which do 
not require correction, as discussed above. Similarly, for the inclusive data (above 2 GeV), the 
resulting VP corrections would be larger without the important recent data from BES which 
are more accurate than earlier measurements and have been corrected appropriately already. 

To estimate the uncertainties in the treatment of VP corrections, we take half of the shifts 
for all channels summed in quadrature^. The total error due to VP is then given by 



/all channels i 

^^vp,excl+incl _ ^ | ^ (^aJ^'Y 1 = ^'^O >< ^O"^" " (18) 




Alternatively, we may assume these systematic uncertainties are highly correlated and prefer 
to add the shifts linearly. For this results in a much smaller error due to cancellations of 

*For data sets with no correction applied, the shifts Aa^P are obviously zero. To be consistent and conser- 
vative for these sets (CLEO, LENA and Crystal Ball) we assign vacuum polarization corrections, but just for 
the error estimate. This results in a total shift of the inclusive data of Aa^P"'"'^'- = —0.94 x 10^^°, rather than 
the -(0.54 + 0.07) x 10"i° implied by TableEl 



9 



the VP corrections, and we prefer to take the more conservative result ()18|) as our estimate of 
the additional uncertainty. However, for Aohad? no significant cancellations are found to take 
place between channels, so adding the shifts linearly gives the bigger effect. Hence for Aohad 
we estimate the error from VP as 

1 all channels i 

dAaZr'^''''' =2 E A(A«had)^P'^ = 1.07 x 10'^ . (19) 
2.2 Final state radiative corrections 

For all the e~^e^ tt+tt^ data (except CMD-2 iT, whose values for c^^^^-, already contain 
final state photons) and e^e~ —>■ K^K~ data, we correct for the final state radiation effects by 
using the theoretical formula 

Cfsr = l + r7(s)a/7r, (20) 

where ri{s) is given e.g. in |l()2p . In the expression for //(s), we take m = m^r for tt+tt^, 
and m = tjik for K^K~ production. Although the formula assumes point-like charged scalar 
bosons, the effects of vr and K structure are expected to be small at energies not too far 
away from the threshold, where the cross section is significant. The above factor corrects the 
experimental data for the photon radiation effects, including both real emissions and virtual 
photon effects. Because there is not sufficient information available as to how the various sets 
of experimental data are corrected for final state photon radiative effects, we include 50% of 
the correction factor with a 50% error. That is, we take 

Q.,= ('l + 0.5„W^^)±0,5„W^^. (21) 

SO that the entire range, from omitting to including the correction, is spanned. The esti- 
mated additional uncertainties from final state photon radiation in these two channels are 
then numerically Sa^^''''^''' = 0.68 x 10"^° and 6a^^''^^^~ = 0.42 x 10"^°, and for Aohad, 
SAa^!^^^^ = 0.04 X 10""^ and 6Aa^^l^^^ = 0.06 x 10""^. For all other exclusive modes we 
do not apply final state radiative corrections, but assign an additional 1% error to the contri- 
butions of these channels in our estimate of the uncertainty from radiative corrections. This 
means that we effectively take 

Ckr = 1 ± 0.01 (22) 
for the other exclusive modes such as Svr, 7r°7, rj'j, An, bn, KKnn, etc., which gives 

^^fer, other _ q ^ (23) 

5A4l'd°'''" = 0.10 X 10"^ (24) 

^For the tt+tt" contribution very close to threshold, which is computed in chiral perturbation theory, we 
apply the exponentiated correction formula (47) of |1U2| . For a detailed discussion of FSR related uncertainties 
in TT+TT" production see also |103| . 
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2.3 Radiative corrections for the narrow {J /ip^ip' ^T) resonances 



The narrow resonance contributions to the dispersion integral are proportional to the leptonic 
widths T{y —>■ e^e~). The leptonic widths tabulated in |104j contain photon vacuum polariza- 
tion corrections, as well as final state photon emission corrections. We remove those corrections 
to obtain the bare leptonic width 

Vl = a,,T{V e+e-) (25) 

where 

l + (3/4)a/7r ^ ^ 

Since a reliable evaluation of a{my) for the very narrow J/ip,ip' and T resonances is not 
available, we use a(— my) in the place of a{my) in ^I^. The correction factors obtained in 
this way are small, namely C^es = 0.95 for J/ip and ip' , and 0.93 for T resonances, in agreement 
with the estimate given in |lU5j . A more precise evaluation of the correction factor will be 
discussed elsewhere |lU6j . 

To estimate the uncertainty in the treatment of VP corrections, we take half of the errors 
summed linearly over all the narrow resonances. In this way we found 

^ y=J/V','(/''.Tr 

= (0.15 + 0.04 + 0.00) X 10"^° (27) 
= 0.19 X 10"^°, (28) 



V=J/tp,tp',T 

4 



= (0.17 + 0.06 + 0.02 + 0.00) x 10"^ (29) 
= 0.25 X 10-^ (30) 

where the three numbers in (j27|) mean the contributions from J /if), if)' and T(15' — 65"), re- 
spectively. Similarly, the four numbers in ()29j) are the contributions from J/ipjip' ,T{1S) and 
T{2S-6S). 



2.4 Combining data sets 

To evaluate dispersion integrals (jD) and (0) and their uncertainties, we need to input the 
function R{s) and its error. It is clearly desirable to make as few theoretical assumptions as 
possible on the shape and the normalization of R{s). Two typical such assumptions are the use 
of Breit-Wigner shapes for resonance contributions and the use of perturbative QCD predictions 
in certain domains of s. If we adopt these theoretical parameterizations of R{s), then it becomes 
difficult to estimate the error of the integral. Therefore, we do not make any assumptions on the 
shape of R{s), and use the trapezoidal rule for performing the integral up to i/i = 11.09 GeV, 
beyond which we use the most recent perturbative QCD estimates, including the complete quark 
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mass corrections up to order a|, see e.g. jlU7j . This approach has been made possible because 
of the recent, much more precise, data on 2tt, Svr, KK, n^'y, rj'y channels in the u and resonant 
regions^°. Although this procedure is free from theoretical prejudice, we still have to address 
the problem of combining data from different experiments (for the same hadronic channel), 
each with their individual uncertainties. If we would perform the dispersion integrals (^, (0) 
for each data set from each experiment separately and then average the resulting contributions 
to (or Aahad), this, in general, would lead to a loss of information resulting in unrealistic 
error estimates (as discussed e.g. in |lU8j ). and is, in addition, impracticable in the case of data 
sets with very few points. On the other hand, a strict point-to-point integration over all data 
points from different experiments in a given channel would clearly lead to an overestimate of the 
uncertainty because the weighting of precise data would be heavily suppressed by nearby data 
points of lower quality. The asymmetry of fluctuations in poorly measured multi-particle final 
states and in energy regions close to the thresholds could in addition lead to an overestimate 
of the mean values of and Aahad- 

For these reasons, data should be combined before the integration is performed. As different 
experiments give data points in different energy bins, obviously some kind of 're-binning' has 
to be applied^^. The bin-size of the combined data will depend, of course, on the available 
data and has to be much smaller in resonance as compared to continuum regimes (see below). 
For the determination of the mean R value, within a bin, the R measurements from different 
experiments should contribute according to their weight. 

The problem that the weight of accurate, but sparse, data may become lower than inaccu- 
rate, but densely-populated, data is well illustrated by the toy example shown in Fig. |21 The 
plots show two hypothetical sets of R data. The set shown by circles has many data points 
with large statistical and a 30% systematic error. The second set has only two data points, 
shown by squares, but has small statistical and only a 1% systematic error. (The length of the 
error bars of each point is given by the statistical and systematic error added in quadrature, 
whereas the little horizontal line inside the bar indicates the size of the statistical error alone.) 
Two alternative ways of treating the data are shown in Fig. |2l together with the respective 
contribution to a^, which follows from the trapezoidal integration. In the first plot, the impact 
of the two accurate data points is local (with a 5 MeV cluster size no combination with the the 
other set takes place and only two of the less accurate points around 1.7 GeV are combined), 
and we see that the integral has a 30% error. In the second plot, we have assumed that R{s) 
does not change much in a 50 MeV interval, and hence have combined data points which lie in 
50 MeV 'clusters'. In this clustering process, the overall normalization factors of the two data 
sets are allowed to vary within their uncertainties. In the toy example, this means that in the 
upper plot no renormalization adjustment takes place, as there is no cluster with points from 
both data sets. In the lower plot, however, the points of the more accurate set 2 are binned 
together in the clusters with mean energies 1.51 and 1.83 GeV and lead to a renormalization of 

^"The J/ip^ip' and the T resonances are still treated in the zero-width approximation. 

^^Another possibility to 'combine' data, is to fit them simultaneously to a function with enough free param- 
eters, typically polynomials and Breit-Wigner shapes for continuum and resonance contributions, see e.g. |99| . 
We decided to avoid any such prejudices about the shape of R and possible problems of separating continuum 
and resonance contributions. 
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all the points of the less accurate set by a factor 1/1.35. (Vice versa the adjustment of set 2 is 
marginal, only (1/0.9995), due to its small errors.) It is through this renormalization procedure 
that the sparse, but very accurate, data can affect the integral. As a result, in the example 
shown, the value of the integral is reduced by about 30% and the error is reduced from 30% 
to 15%. The goodness of the fit can be judged by the Xmm P^r degree of freedom, which is 
0.61 in this toy example. We find that by increasing the cluster size, that is by strengthening 
our theoretical assumption about the piecewise constant nature of R, the error of the integral 
decreases (and the xLin P^r degree of freedom rises). Note that the 'pull down' of the mean 
R values observed in our toy example is not an artefact of the statistical treatment (see the 
remark below) but a property of the data. 

More precisely, to combine all data points for the same channel which fall in suitably chosen 
(narrow) energy bins, we determine the mean R values and their errors for all clusters by 
minimising the non-linear function 



x'iRm, fk) = E [(1 - h) l^hf + E E [(^f '"^^ - hRra) I ^R 



{fc,m} 



k=l 



m=l 



(31) 



Here Rm and fk are the fit parameters for the mean R value of the rri^^ cluster and the overall 
normalization factor of the k^^ experiment, respectively. r\^'^^ and d-R^'^''""'^ are the R values 
and errors from experiment k contributing to cluster m. For d-Rl^'™^ the statistical and, if given, 
point-to-point systematic errors are added in quadrature, whereas d/^ is the overall systematic 
error of the k^^ experiment. Minimization of (|31|) with respect to the (A''exp + -/Vciust) parameters, 
fk and Rmi gives our best estimates for these parameters together with their error correlations. 

In order to parameterize R{s) in terms of Rm-, we need a prescription to determine the 
location of the cluster, ^/s = Em- We proceed as follows. When the original data points, which 
contribute to the cluster m, give 



R{V^ = K 



{fc,r7i}\ 



{k,m} 



+ {dfk 



from the kth experiment, we calculate the cluster energy by 



Er, 



{k,m.} 



E E 



1 



(di?f'"^>) +(d/.)^ 



-E„ 



{k,m,} 



N, 



{ fc ,m} 



E E 



(di?f'"^^) +(d/. 



(32) 



(33) 



where the sum over k is for those experiments whose data points contribute to the cluster m. 
Here we use the point-to-point errors, dRj'''"^\ added in quadrature with the systematic error, 
dfk, to weight the contribution of each data point to the cluster energy Em- Alternatively, we 
could use just the statistical errors to determine the cluster energies Em- We have checked that 
the results are only affected very slightly by this change for our chosen values for the cluster 
sizes. 

The minimization of the non-linear function with respect to the free parameters Rm and 
fk is performed numerically in an iterative procedure^^ and we obtain the following parameter- 
ization of R{s)-- 

R{s = E l) = Rm = Rm±dRm, (34) 

^^Our non- linear definition H31|l of the function avoids the pitfalls of simpler definitions without rescaling 
of the errors which would allow for a linearized solution of the minimization problem, see e.g. |1(J9I IllU] . 
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Figure 2: Two toy data sets chosen to illustrate the problems of combining precise with less 
precise data. The upper plot shows the result obtained with a very small 'cluster' size. The 
lower shows the data clustered in 50 MeV bins, which allows renormalization of the data within 
their systematic errors. Here the (much less precise) points of set 1 are renormalized by 1/1.35 
whereas the two precise points of set two are nearly unchanged (1/0.9995). The length of the 
error bars give the statistical plus systematic errors added in quadrature for each data point. 
The small horizontal lines in the bars indicate the size of the statistical errors. The error band 
of the clustered data is defined through the diagonal elements of the covariance matrix. 
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where the correlation between the errors dRm and di?„, 

Pcorr("^,n) = V{m,n)/{dRm)idR„), (35) 
with V{m,m) = {dRm)^, is obtained from the covariance matrix V{m,n) of the fit, that is 

= xLn + E IliRm - Rm)V~\m, n)(/?„ - Rn). (36) 

m=l n=l 

Here the normahzation uncertainties are integrated out. We keep the fitted values of the 
normalization factors fk 

fk = fk. (37) 

The function takes its minimum value xLin when R^ = Rm and fk = fk- The goodness of 
the fit can be judged from 

A mm A mm ^'QS^ 

d:^ ~ Efc(iVfc-l)-iVcluster' ^ ^ 

where Y^k^k stands for the total number of data points and Y^k{~^) stands for the overall 
normalization uncertainty per experiment. Once a good fit to the function R{s) is obtained, 
we may estimate any integral and its error as follows. Consider the definite integral 

= / dsR{s)K{s) = 2 dEER{E^)K{E^) = I±AI. (39) 
When a = E^ < En = b, the integral / is estimated by the trapezoidal rule to be 

~-3lLEmR^Kn, + ^" 'f'^-' EnRnKn + E ^'^'~'^'~' EkRkKk] ,m 

^ k=m+l ^ I 

where = K{El), and its error A/ is determined, via the covariance matrix V, to be 

(A'f ^ f Ej|mof (41) 

k=m l=m ' 
n 

= ii^k+i - Ek^i)Ek Kk ) Vik, I) ((E,+i - Ei.,)EiKi) , (42) 

k,l=m 

where Em-i = Em and En+i = En at the edges, according to (pH) . When the integration 
boundaries do not match a cluster energy, we use the trapezoidal rule to interpolate between 
the adjacent clusters. 

We have checked that for all hadronic channels we find a stable value and error for a^'^'^'^'-', 
together with a good^^ fit if we vary the minimal cluster size around our chosen default 
values (which are typically about 0.2 MeV for a narrow resonance and about 10 MeV or larger 
for the continuum). For the most important 7i~^7t~ channel we show in Fig. IHlthe behaviour of 

^^However, there are three channels for which Xmin/'^-'^-^- > l-^? indicating that the data sets are mu- 
tually incompatible. These are the e+e" 7r+7r~7r+7r^, tt+tt^tt^, 7r+7r~7r°7r° channels with Xmin/^-'^-^- ~ 
2.00, 1.44, 1.28 respectively. For these cases the error is enlarged by a factor of ^x^^j^^/d.o.f.. Note that for 
the four pion channel a re-analysis from CMD-2 is under way which is expected to bring CMD-2 and SND data 
in much better agreement |111| . If we were to use the same procedure, but now enlarging the errors of the data 
sets with Xmin/d-O.f. > 1, then we find that the experimental error on our determination of Qrhad.LO increased 
by less than 3% from the values given in H125|) and H12t)f) below. 
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the contribution to a^, its error and the quahty of the fit expressed through x^j^/d.o.f. as a 
function of the typical cluster size 6. It is clear that very large values of S, even if they lead to 
a satisfactory Xmin) should be discarded as the fit would impose too much theoretical prejudice 
on the shape of R{s). Thus, in practice, we also have to check how the curve of the clustered 
data, and its errors, describe the data. One would, in general, try to avoid combining together 
too many data points in a single cluster. 

In Table |3] we give the details of the clustering and non-linear fit for the most relevant 
channels. The fits take into account data as cited in Table ^ with energy ranges as indicated 
in the second column of Table |3 We use clustering sizes 6 as displayed in the third column. 
In the TT^Ti^Ti^, K^K~ and K^K^ channels the binning has to be very fine in the uj and 
resonance regimes; the respective values of the clustering sizes in the continuum, {uj and) 
regions are given in the Table. The Xmin/^-o.f. displayed in the fourth column is always good, 
apart from the three channels vr+vr^Tr+Tr^, tt+tt^tt^ and 7r+7r~27r'', in which we infiate the error 
as mentioned above. In most cases the fit quality and result is amazingly stable with respect to 
the choice of the cluster size, indicating that no information is lost through the clustering. Table 
13] also gives information about the contribution of the leading channels to within the given 
ranges. For comparison, the last column shows the contributions to obtained by combining 
data without allowing for renormalization of individual data sets through the fit parameters fk- 
In this case, we use the same binning as in the full clustering, but calculate the mean values 
Rm just as the weighted average of the R data within a cluster: 



(These Rm values are actually used as starting values for our iterative fit procedure.) The point- 
to-point trapezoidal integration ()4()j) with these Rm values from without the fit neglects 
correlations between different energies. As is clear from the comparison of columns six and 
eight of Table |3J such a procedure leads to wrong results, especially in the most important 
TT+vr" channel. 

As explained above, the dispersion integrals (jH) and (0) are evaluated by integrating (us- 
ing the trapezoidal rule ()4U|) for the mean value and (jl^ for the error and thus including 
correlations) over the clustered data directly for all hadronic channels, including the u and 
resonances. Thus we avoid possible problems due to missing or double-counting of non-resonant 
backgrounds. Moreover interference effects are taken into account automatically. As an exam- 
ple we display in Fig. [7| the most important tt+tt^ channel, together with an enlargement of 
the region of p-uo interference. As in Fig. the error band is given by the diagonal elements 
of the covariance matrix of our fit, indicating the uncertainty of the mean values. Data points 
are displayed (here and in the following) after application of radiative corrections. The error 
bars show the statistical and systematic errors added in quadrature and the horizontal markers 
inside the error bars indicate the size of the statistical error alone. 

In the region between 1.43 and ~ 2 GeV we have the choice between summing up the exclu- 
sive channels or relying on the inclusive measurements from the 772, MEA, M3N and ADONE 
experiments |H3]^|HZ1- Two-body final states were not included in these analyses. Therefore we 
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Figure 4: The behaviour of R obtained from inclusive data and from the sum of exclusive 
channels, after clustering and fitting the various data sets. Note the suppressed zero of the 
vertical scale. 

correct the R data from 772, MEA and ADONE for missing contributions from tt+tt", K^K~ 
and K^K^l, estimating them from our exclusive data compilation.^^ The corrections are small 
compared to the large statistical and systematic errors and energy dependent, ranging from 
up to 7% at 1.4 GeV down to about 3% at 2 GeV. In addition, we add some purely neutral 
modes to the inclusive data, see below. Surprisingly, even after having applied these correc- 
tions, the sum of exclusive channels overshoots the inclusive data. The discrepancy is shown 
in Fig. U where we display the results of our clustering algorithm for the inclusive and the 
sum of exclusive data including error bars defined by the diagonal elements of the covariance 
matrices (errors added in quadrature for the exclusive channels). We study the problem of this 
exclusive/inclusive discrepancy in detail in Section |3 



3 Evaluation of the dispersion relations for ajf ' and 

Aahad 

Here we use dispersion relations (jH) and (0) to determine a^^'^'^'~' and Aa;had(^|) respectively^^, 
which in turn we will use to predict g — 2 oi the muon (in Section [7j) and the QED coupling 

^^We do not correct the data from M3N as they quote an extra error of 15% for the missing channels which 
is taken into account in the analysis. 

-^^It is conventional to compute Aahad for 5 quark flavours, and to denote it by Aa[^^^. For simplicity of 
presentation we often omit the superscript (5), but make the notation explicit when we add the contribution of 
the top quark in Section |S1 
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a{M'^) (in Section |HJ. The dispersion relation (0} has the form 

where (TYiadi^) total cross section for e+e^ hadrons (+7) at centre-of-mass energy y^, 

as defined in (j2)). For s > 4m^ the kernel function K{s) is given by |112j 

X(s > 4my = 1^(2-2; ) + ^ ^5 ^ (^ln(l + x) - a; + y j + j—^x^lnxj ,(45) 

with X = (1 — Pfj)/{1 + /3^) where = ^Jl — ^m?^/ s; while for s < ^m?^ the form of the kernel 
can be found in and is used to evaluate the small 7r°7 contribution to aJ^^^'LO Dispersion 

relation (jSJ, evaluated at s = M|, may be written in the form 

A„,..,M|,^-^p£d.^. (46) 

To evaluate (jiij) and (jlUj) we need to input the function cr°a^d(s) and its error. Up to y/s ~ 
2 GeV we can calculate (J^^^ from the sum of the cross sections for all the exclusive channels 
e'^e" — > TT+TT^, 7r+7r^7r°, etc. On the other hand for y/s > 1.4 GeV the value of a'^^ can be 
obtained from inclusive measurements of e^e~ —>■ hadrons. Thus, as mentioned above, there is 
an 'exclusive, inclusive overlap' in the interval 1.4 < ^/s < 2 GeV, which allows a comparison of 
the two methods of determining a^^^^ from the data. As we have seen, the two determinations 
do not agree, see Fig. EJ It is worth noting that the data in this interval come from older 
experiments. The new, higher precision, Novosibirsk data on the exclusive channels terminate 
at s/s ~ 1.4 GeV, and the recent inclusive BES data [HHllHn] start only at a/s ~ 2 GeV. Thus 
in TableElwe show the contributions of the individual channels to qJ^'^'^'LO g^j^^ Aahad(Af|) using 
first inclusive data in the interval 1.43 < y/s < 2 GeV, and then replacing them by the sum of 
the exclusive channels. 

Below we describe, in turn, how the contributions of each channel have been evaluated. 
First we note that narrow u and contributions to the appropriate channels are obtained by 
integrating over the (clustered) data using the trapezoidal rule. We investigated the use of 
parametric Breit-Wigner forms by fitting to the data over various mass ranges. We found no 
significant change in the contributions if the resonant parameterization was used in the region 
of the u and peaks, but that the contributions of the resonance tails depend a little on the 
parametric form used. The problem did not originate from a bias due to the use of the linear 
trapezoidal rule in a region where the resonant form was concave, but rather was due to the 
fact that different resonant forms fitted better to different points in the tails. For this reason 
we believe that it is more reliable to rely entirely on the data, which are now quite precise in 
the resonance regions. 

3.1 7r^7 channel 

The contribution of the e+e^ —>■ tc^'j channel defines the lower limit, i/sth ~ ^-^^ 
dispersion integrals. There exist two data sets pT| |S2j for this channel, which cover the interval 
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0.60 < ^/s < 1.03 GeV (see Fig. El). After clustering, a trapezoidal rule integration over this 
7r°7 energy interval gives a contribution 

a^(7r°7, 0.6 < < 1.03 GeV) = (4.50 ± 0.15) x 10"^° (47) 

and 

Aahad(7r°7, 0.6 < < 1.03 GeV) = (0.36 ± 0.01) x 10"^ (48) 

In Fig. El we show an overall picture of the e+e^ 7r°7 data and a blow up around the p-u 
region. 

The use of the trapezoidal rule for the interval m,r < \^ < 0.6 GeV would overestimate 
the contribution, since the cross section is not linear in ^/s. In this region we use chiral 
perturbation theory (ChPT), based on the Wess-Zumino-Witten (WZW) local interaction for 
the 7r''77 vertex, 

Cwzw = -^vr°e'^^"'^F^,FA., (49) 
with — 93 MeV, which yields 

( + - o^ _ 8mrr(7r0_^27) / m^' 

Since the electromagnetic current couples to 7r°7 via u meson exchange, the low-energy cross 
section can be improved by assuming the w-meson dominance |113j . which gives 

avMD(e+e- ^ 7r%) = apt(e+e- ^ 7r%) (^^] • (51) 

s J 

We find 

S(^°7' Vs < 0.6 GeV) = (0.13 ± 0.01) x 10"^°, (52) 

while the contribution to Aohad is less than 10~^. The agreement of the prediction of 1)511) for 
the 7r°7 cross section with the SND data just above 0.6 GeV is shown in Fig. IHl 



3.2 TT+TT channel 

We use 16 data sets for e^e~ tt+tt" jTHI, [lEj^lSD! which cover the energy range 0.32 < ^/s < 
3.0 GeV. Some older data with very large errors are omitted. In Fig. [3 we show the region 
around p, which gives the most important contribution to (7 — 2 of the muon. 

The 7r+7r~ contributions^^ to a^^^d.LO ^-^^ Aahad(^|)5 obtained by integrating clustered data 
over various energy intervals, are shown in Table IHl As seen from the Table, if we integrate 

^^If we were to leave out the dominant tt+tt" data from CMD-2 altogether, we find 491.33 ± 8.47, instead 
of 503.38 ± 5.02, for the Tr+Tr" contribution from the interval 0.32 < < 2 GeV. (The Xmi„/d-0-f- of the fit 
which clusters the data would be even slightly better, 1.00 instead of 1.07.) This means that after re-analysis 
the CMD-2 data dominate the error but do not pull down the contribution, but rather push it up! 
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Figure 5: Data for a^{e~^e —>■ 7r°7). The shaded band shows the behaviour of the cross section 
after clustering and fitting the data. The second plot is an enlargement in the region of the oo 



resonance. 
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Figure 6: Predictions for a^{e~^e~ ir^'y) from ChPT compared with low energy experimental 
data from the SND collaboration [31]. In the figure we have multiplied m'j^/{3s)K{s) by the 
cross section so that the area below the data is proportional to the contribution to a^. The 
continuous curve, which is obtained assuming Vector Meson (u) Dominance (VMD), is used 
for s < 0.36 GeVl 
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Figure 7: e^e^ tt+tt^ data up to 1.2 GeV, after radiative corrections, where the shaded band 
shows the result, cr''(7r+7r~) (obtained from Rm of (jSH)), of our fit after clustering. The width 
of the band indicates the error on the a^{7T~^7T~) values, obtained from the diagonal elements 
of the full covariance matrix. The second plot is an enlargement of the p-u interference region. 
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over the data up to 1.43 GeV, we obtain 

a^{n+n-, 0.32 < < 1.43 GeV) = (502.78 ± 5.02) x 10'^°, (53) 

Aahad(7r+7r~, 0.32 < < 1.43 GeV) = (34.39 ± 0.29) x 10"^ (54) 

If we integrate up to 2 GeV, instead of 1.43 GeV, we obtain 

a^(7r+7r", 0.32 < < 2 GeV) = (503.38 ± 5.02) x 10"^°, (55) 

Aahad(vr+7r~, 0.32 < < 2 GeV) = (34.59 ± 0.29) x 10"^ (56) 



The contribution of the vr+vr" channel is dominated by the p-meson, and hence the differences 
between Eqs. (|^ and (|^ is smalL If we use the CMD-2 data before the recent re-analysis pUj . 
we have 

a^in+n-, 0.32 < < 1.43 GeV, old CMD-2 data) = (492.66 ± 4.93) x 10"^°, (57) 
Aahad(vr+7r-, 0.32 < < 1.43 GeV, old CMD-2 data) = (33.65 ± 0.28) x 10"^ (58) 

The comparison of ()53|) and ()57|) shows the effect of the re-analysis of the recent CMD-2 data, 
which is an upward shift of the central value by roughly 2% in this interval. 

It is interesting to quantify the prominent role of these most precise CMD-2 e+e^ tt+tt" 
data, which have a systematic error of only 0.6%. If we were to omit these CMD-2 data in the 
central p regime altogether, the contribution of this channel to would decrease by roughly 
12.1 X 10~^°, i.e., by ~ 2.4%, whereas the error would increase by about 3.4 x 10~^°, i.e., by 
~ 68% in the interval 0.32 < ^/s < 1.43 GeV. 

In the threshold region, below 0.32 GeV, we use chiral perturbation theory, due to the lack 
of TT+vr" experimental data. The pion form factor Ft^{s) is written as 

F^{s) = 1 + ]:{r^)^s + c^s' + 0{s'), (59) 
6 

with coefficients determined to be |114j 

(r^)^ = 0.431 ± 0.026 (fm^), = 3.2 ± 1.0 (GeV"^), (60) 

by fitting to space-like pion scattering data |115j . Fig. |H1 compares the prediction with the 
(time-like) experimental data which exist for ^/s > 0.32 GeV. The contributions from the 
threshold region are 

a^(7r+7r-,v^< 0.32 GeV) = (2.36 ± 0.05) x 10"^°, (61) 
Aahad(7r+7r", < 0.32 GeV) = (0.04 ± 0.00) x 10'^, (62) 

and are also listed in the last row of Table IHl Though these contributions are small, for it is 
non- negligible. 

In the calculation of the contribution from the threshold region, we have included the effect 
from final state (FS) radiative corrections. In Ref. (102j both the 0{a) correction and the 
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Figure 8: The e~^e~ tt+tt^ data near the threshold. In the figure we have multiphed 
m'^/{3s)K{s) by the cross section so that the area below the data is proportional to the con- 
tribution to a^. The theoretical curve obtained from chiral perturbation theory is also shown, 
and is used up to = 0.32 GeV {s = 0.10 GeV^). 

exponentiated formula for the FS corrections are given. If we do not apply the FS correction, 
we would obtain 

a^(7r+7r", < 0.32 GeV) = (2.30 ± 0.05) x 10"^°. (63) 

However, if we include the FS corrections, we have 

a^{n+n-, ^ < 0.32 GeV, 0{a) FS corr.) = (2.36 ± 0.05) x 10"^°. (64) 

We obtain the same contribution if we use the exponentiated formula, which we have used in 
all the tables in the paper. The effect of final state radiation is to increase the contribution by 
about 3 %, whether the 0{a) or the exponentiated form is used. Similarly, the contribution 
from this region to Aahad(^I) is given by 

Aahad(7r+7r", < 0.32 GeV, exponentiated FS corr.) = (0.04 ± 0.00) x 10-^ (65) 

so here the contribution from the threshold region is totally negligible. 

3.3 TT+TT TT*^ channel 

We use ten experimental data sets for the tt+tt^tt^ channel PHl El 1221 01], |HZ]~|12|, which 
extend up to 2.4 GeV, although the earlier experiments have large errors, see Fig. 1^1 Since 
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the data for this channel are not very good, we inflate the error by a factor of -y/xmin/'i-o.f., 

which is 1.20 for this channel. (We inflate the error by a factor of ^Xmm/d-o.f. whenever 
Xmin/d-O-f- > 1-2, as dlscusscd in Section see Table El) We discard the data points below 
0.66 GeV, in favour of the predictions of chiral perturbation theory jll6tlll7] . see Fig. ^| The 
contributions to a^'^'^'^'~' and Aahad(^l) are 

a^(7r+7r-7r°, 0.66 GeV < ^< 1.43 GeV, data) = (46.43 ± 0.90) x 10"^°, (66) 
Aahad(vr+7r-7r°, 0.66 GeV < < 1.43 GeV, data) = (4.33 ± 0.08) x 10^^ (67) 

respectively. 

In the threshold region, below 0.66 GeV, we use chiral perturbation theory |116| IllTj . due 
to the lack of good tc^tt'tt^ experimental data, see Figs. Eland^J The contributions to a^^'^'^'^ 
and Aahad(^I) froKi the threshold region are 

a^(7r+7r"7r°, < 0.66 GeV, ChPT) = (0.01 ± 0.00) x 10~^°, (68) 
Aahad(vr+7r"7r°,v^< 0.66 GeV, ChPT) = (0.00 ± 0.00) x 10"^ (69) 

There is a tendency that the ChPT prediction with the tu-dominance undershoots the lowest- 
energy data points. Because of the smallness of the threshold contribution, we do not attempt 
further improvement of the analysis. 

3.4 ?77 channel 

We use five data sets from SND |S21ESI and CMD-2 [311123 EE]- We divide the data set given 
in Ref. jHEl into two parts at 0.95 GeV since it has different systematic errors below and above 
this energy. 

Since the lowest data point starts only at 690 MeV, we use ChPT at the threshold region up 
to the lowest-energy data point. We summarize our method in Appendix A, according to which 
the contribution from the region to a^^'^'^'^ is less than 10~^^, which can be safely neglected. 
The contribution to Aohad is also small, and less than 10"''. In Fig. ^2 we show the threshold 
region of the 777 production cross section and our prediction from ChPT. 

Above the lowest-energy data point we integrate over the data. In Fig.^Jwe show the overall 
picture of the 777 production cross section and our result of the clustering. After integrating 
over 0.69 < y/s < 1.43 GeV we obtain 

0^(777, 0.69 < < 1.43 GeV) = (0.73 ± 0.03) x 10~^°, (70) 
Aahad(^7, 0.69 < ^/s < 1.43 GeV) = (0.09 ± 0.00) x 10"^ (71) 

3.5 47r, Stt, Gtt and tjtt^tt^ channels 

For the An channel, we have data for the 27r+27r^ and the 7r+7r^27r° final states. (The reaction 
e~^e~ — > 7* 47r° is forbidden from charge conjugation symmetry.) 
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Figure 9: The data for (T°(e"'"e~ ti^ti^tt^) together with an expanded version in the u and 
resonance regions. The shaded band shows the result of our fit after clustering. In the analysis 
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For the 27r+27r- channel, we use thirteen data sets [22l IHOl EH , [HH|-|nn], [H2]-|M!, see Fig. 

IT^ Since the data for this channel are not very consistent to each other, we inflate the error 
by a factor of 



A mm/ 



yd.o.f. = 1.41. We note, in particular, that the compatibility between the 
data from SND and CMD-2 is poor. This may improve after the re-analysis of the CMD-2 data 
for this channel is completed jlllj . The contribution from this channel is 



a^(27r+27r", ^ < 1.43 GeV, data) 



(6.16 ±0.32) X 10"^°, 
(1.27 ±0.07) X 10"^ 



Aahad(27r+27r", < 1.43 GeV, data) = 
For the 7r+7r^27r'' channel, we use eight data sets |Sn]-[HZl, see Fig. ^[ which contribute 



(72) 
(73) 



a^(7r+7r"27r°, < 1.43 GeV, data) 
Aahad(7r+7r-27r°, < 1.43 GeV, data) 



(9.71 ±0.63) X 10"^°, 
(1.86 ±0.12) X 10"^ 



(74) 
(75) 



A mm 



/d.o.f. = 1.13 as discussed in 



For the tt+tt 27r° channel we have inflated the error by 
Section in 

For the Svr channel, there exist data for the 27r+27r^7r° and the 7r+7r^37r° final states. (The 
reaction e^e~ ^ 7* — > Stt^ is forbidden from charge conjugation symmetry.) We use five data 
sets for the 27t^27i^tt^ channel |221 EOl ESI EZl E21 , and one data set for the 7r+7r^37r'' channel |50j . 
We integrate over the clustered data, which gives 



and 



a^(27r+27r-7r°, < 1.43 GeV, data) 
Aahad(27r+27r"7r°, < 1-43 GeV, data) 



a/,(7r+7r"37r°,v^< 1.43 GeV, data) 
Aahad(7r+7r"37r°,A/s < 1.43 GeV, data) 



(0.26 ±0.04) X 10~^°, 
(0.06 ±0.01) X 10-^ 



(0.09 ±0.09) X 10"^°, 
(0.02 ±0.02) X 10-^ 



(76) 
(77) 



(78) 
(79) 
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respectively. For the Svr channels we do not inflate the error since the Xmm/'^-O-^- values are 

xLn/d.o.f.(27r+27r-7rO) = 0.90, (80) 
xL„/ci.o.f.(7r+7r-37rO) = 1.07. (81) 

For the Gn channel, there are data for the Stt+Stt" and the 27r+27r^27r'^ flnal states, but not 
for the tc^ti^Ati^ final state. For the tt^ti^Att^ channel we estimate the contribution to and 
^Q^had by using an isospin relation. The reaction e+e" — > 7* ^ 67r° is forbidden from charge 
conjugation. 

We use four data sets for the Stt+Stt" channel [501 IZH IZ21 • M3N [SUl provides the lowest 
data point which is at 1.35 GeV, which we do not use since it has unnaturally large cross 
section with a large error, (1.56 ± 1.11) nb, compared with the next data point from the same 
experiment, (0.10 ± 0.31) nb at 1.45 GeV. The first data points from CMD and DM1 [H] 
contain data with vanishing cross section with a finite error, which result in points with zero 
cross section even after clustering. We do not use such points when integrating over the data. 
Thus the first data point after clustering is at 1.45 GeV. Our evaluation of the contribution 
from the 37r"*"37r~ channel from the region ^/s < 1.43 GeV is zero for both and Aa^ad- 

For the 27r+27r^27r° channel we use five data sets [301 13^ I37| 1^ I72j . which cover the energy 
interval from 1.32 GeV to 2.24 GeV. The trapezoidal integration gives us 

a^(27r+27r"27r°,v^< 1.43 GeV, data) = (0.12 ± 0.03) x 10"^°, (82) 
Aahad(27r+27r"27r°,v^< 1.43 GeV, data) = (0.03 ± 0.01) x 10"^ (83) 

For the tt^tt'^Atc^ channel we use the multipion isospin decompositions |118t I119j of both 
the e~^e~ Gvr channel and the r GirUr decays, which are summarised in the Appendix of 
Ref. |120j . Then using the measured ratio |121j of r~ —>■ 27r~7r+37r°z/^ and r~ —>■ 37r^27r+7r°z/T- 
decays, and the observed u dominance of final states of r ^ GnUr decays [1203, 

a(7r+7r"47r°) = 0.031 (T(27r+27r"27r°) + 0.093 (T(37r+37r"). (84) 

Hence we obtain the small n'^ir'Air^ contribution^^ shown in Table We assign a 100 % error 
to the cross section computed in this way. For a}^^'^'^^ and Aohad they are less than 10~^^ and 
10~^, respectively, when integrated up to 1.43 GeV. 

For the 7]7r~^7T~ channel, we use two data sets [011 IZS]- The entry for the rjix^Ti^ channel in 
Table 101 shows the contribution of cr(e^e~ ri7T~^7r~) multiplied by (1 — -B(?7 37r°) — B{ri — >■ 
TT^n^Ti^)) ~ 0.448, since these rj decay modes are already included in the contribution of the 
57r channels. The contributions to the muon g — 2 and Aa^ad are 

a^ivi^ 7r°7)7r+7r-, < 1.43 GeV) = (0.07 ± 0.01) x 10"^°, (85) 
Aahad(^(-^ 7r°7)7r+7r~, < 1.43 GeV) = (0.02 ± 0.00) x 10"^ (86) 

"'^^Relation (|84|l was not used in our previous analysis 12 . As a consequence, the (weaker) isospin bound 
then gave a larger contribution for the 7r"''7r~47r*' channel. However DEHZ ^ did use the observed information 
of r — > Girvr decays to tighten the isospin bound. 
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3.6 K^K and KgKi contributions 

For the K^K' channel, we use ten data sets [2112312111211, |13-|1Z1, which extend from 1.0 
GeV to 2.1 GeV, see Fig. HH When integrated, this channel contributes to the muon g — 2 and 
Aahad an amount 

a^(is:+fs:",v^< 1.43 GeV, data) = (21.62 ± 0.76) x 10"^°, (87) 
Aahad(i^+^~,v^< 1-43 GeV, data) = (3.01 ± 0.11) x 10"^ (88) 

For the K^K^ channel, we use ten data sets [101 121 113 HH HI] , which also extend from 
1.0 GeV to 2.1 GeV, see Fig. ^1 Using the trapezoidal rule, the channel gives a contribution 
to the muon g — 2 and Aahad of 

a^(is:^is:^, < 1.43 GeV, data) = (13.16 ± 0.31) x 10"^°, (89) 
Aahad(i^s^°, ^ < 1-43 GeV, data) = (1.76 ± 0.04) x 10"^ (90) 

This channel is the one case where the use of the trapezoidal rule may overestimate the resonance 
contribution, due to the lack of data in certain regions of the resonance tails, see Fig. ^1 
We find that the use of a smooth resonance form in the tails decreases the contributions to 
and Aohad by about 0.15 x 10^^° and 0.02 x 10~"^ respectively. We have therefore increased the 
error in (j89|) and (|90|) to include this additional uncertainty. 

3.7 if A' + nTT contributions 

We take into account the KK + nvr final states for n = 1 and 2. 

For the KKti, in addition to the data for the K^ii^K^ Hi [73 HH] and K+R-n^ [H [7S] 
channels, we use the equalities a{K^7TK) = a{Kg7rK) and a{KgK^7T^) = a{K^K~7r^), which 
follow directly from isospin. The contribution from the K^n^K^ + K^tt^K^ channel is 

a^{K°s7r^K^ + KItt^K^, < 1.43 GeV, data and isospin) = (0.10 ± 0.04) x 10"^°, (91) 
Aahad(i^5vr^ir^ + Kln^K^, < 1.43 GeV, data and isospin) = (0.02 ± 0.00) x 10"^ (92) 

For the ti^ + K'^K^n^ channel, the contribution from the region ^/s < 1.43 GeV is taken 

to be zero since the first data point is at 1.44 GeV. 

To evaluate the KKtttt contribution we use the inclusive data for Ks^ [ZZl, together with 
the cross section relation 

2KsX = KsX + KlX 

= 2KsKl + 2{KsKl + KsKs + KlKl) (tt + tttt) 

+ {Ks + Kl) {Ktx + Ktitx), (93) 

where 2KsX stands for 2ij^{e^e~ KsX) and similarly for the other abbreviations. On the 
right-hand-side vrvr stands for tt^tt^ or tt^tt^. Kit for K^n^ or K~it~^, and Kmr for K~^7i:~7t^ or 
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in certain regions of the (p resonance tails, see the discussion in the text. 



34 



K 7r+7r°. On the other hand, the KKtttt cross section is given by 

KKrni = {KsKl + KsKs + KLKL)in7i) + {Kg + KL){Kn7i) + {K+K-){n7i) 
= 2KsX - 2KsKl - {KsKl + KsKs + KlKl) (27r + tttt) 

-2Ks{Kt() + {K+K'){'n'K) 
= 2{KsX - KsKl- K+K-n - Ks{K7t)), (94) 

where to obtain the second equahty we have used ()93|) . In other words, the total KKtttt contri- 
bution is obtained from twice the inclusive KsX cross section by subtracting the appropriate 
KK and KKtt contributions. For this channel, the contribution from the region ^/s < 1.43 
GeV is also taken to be zero since the data of the K^X final state start from 1.44 GeV. 

3.8 Unaccounted modes 

We still have to take into account contributions from the reactions e+e^ —>■ ujtt^ and e~^e~ — *■ 
cuTT+TT^, in which the u decays radiatively into We used seven data sets for the e~^e — > ujtc 
channel [221 El El, ISHl-EIl, and three data sets jSHl EOI for the e+e" lutt+tt- channel. 
Note that the contributions from the u!{^ 7r+7r^7r°)7r° and uj{^ 71^71^)71^ channels are already 
included as a part of the multi-pion channels. We therefore need simply to multiply the original 
cross section cr(e^e~ un^) by the branching ratio B{uj 7i^'~f) = 0.087 |104j . The same 
comments apply for the un^TT' channel. The two channels give contributions 

a^{uj{^ 7r°7)7r°, < 1.43 GeV) = (0.64 ± 0.02) x 10~^°, (95) 
Acthad(cc'(-^ 7r%)7r°, < 1.43 GeV) = (0.12 ± 0.00) x 10'^, (96) 

and 

a^(cj(^ 7r°7)7r+7r-, < 1.43 GeV) = (0.01 ± 0.00) x 10'^°, (97) 
Aahad(^(-^ 7r%)7r+7r-, < 1.43 GeV) = (0.00 ± 0.00) x 10-^ (98) 

respectively. 

Purely neutral contributions from the direct decays of p and u to 7r°7r°7 can be safely 
neglected, as the branching fractions are of the order 5 x 10~^ and 7 x 10~^ respectively |104[ 
I122| I123j . and are suppressed compared to the decays into 7r°7. 

For the resonance we have so far accounted for the (p ~^ K~^K~, KgK^, Stt, 777 and 
7r°7 channels. Since the branching fractions of these final states add up to 99.8% |104j . we 
must allow for the 0.2% from the remaining final states. To do this, we first note that the 
contribution to a^^'^'^'^ from the K^K~ channel in the region is 

a^(0 ^ K^K~; 2m^+ < < 1.03 GeV) = 16.15 x 10"^°. (99) 

Using this, we estimate that the total contribution from the (p to be 

a^(0) = a^(0 ^ K+K-)/B{(j) K+K') = 32 x 10"^°. 
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Hence we include the small residual contribution 

a^(0 remaining channels) = a^(0) x 0.002 = 0.06 x 10~^°, (100) 

and assign to it a 100% error. In a similar way the contribution AQ;had(0 K'^K^) = 
2.12 X 10~^ is used to estimate 

AQ;had(0 — ^ remaining channels) = 0.01 x 10""^, (101) 

to which we again assign a 100% error. 

3.9 Baryon-pair contribution 

If we are to integrate up to high enough energy to pair-produce baryons, we have to take into 
account the pp and nn final states. The data come from the FENICE [78, 79J, DM1 [H2] and 
DM2 [HUl IHT] collaborations for the pp channel, and from the FENICE collaboration [7S1 US] 
for the nn channel. They do not contribute when we integrate over the exclusive channels only 
up to 1.43 GeV, but if we integrate up to 2.0 GeV, the pp channel gives a contribution of 



a/.(m Vs < 2.0 GeV) = (0.04 ± 0.01) x 10"^°, (102) 

Aahad(pp, Vs < 2.0 GeV) = (0.02 ± 0.00) x 10~^ (103) 
while the nn channel gives 

a^{nn, < 2.0 GeV) = (0.07 ± 0.02) x 10"^°, (104) 

Aahad(^n, < 2.0 GeV) = (0.03 ± 0.01) x 10^^ (105) 



3.10 Narrow resonance {J/i{j,i/j\T) contributions 

We add the contributions from the narrow resonances, J/ip^ip' and T{1S — QS). We treat them 
in the zero-width approximation, in which the total production cross section of a vector meson 
V {V = J/^,^',T) is 

a(e+e- ^V) = l2Ti^^5{s - M^). (106) 

Here F^g is the bare leptonic width of V , 

F", = ae.F(y ^ e+e-), (107) 

where 

{a/a{ml)f 
1 + (3/4)a/7r 

which is about 0.95 for J /if) and ip', and about 0.93 for the six T resonances |105j . We use the 
values compiled in RPP for the leptonic widths, T{y e^e"), and obtain the contributions 

a^{J/i)) = (5.89 ± 0.41) X 10"^°, (109) 

a^{iP') = (1.41 ± 0.12) X 10"^°, (110) 

a^{T{lS)) = (0.05 ± 0.00) X 10"^°, (111) 

a^(T(25-65)) = (0.05 ± 0.00) X 10-^°, (112) 
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Figure 16: Data for the measurement of the inclusive hadronic cross section below 2 GeV (left) 
and above 2 GeV (right). The shaded band shows the behaviour of the hadronic R ratio after 
clustering and fitting the data. 



and 



Aahad(^/^) = (6.65 ± 0.47) x 10^^ (113) 

Aahad(V'') = (2.25 ± 0.19) X 10-^ (114) 

Aahad(T(l^)) = (0.54 ± 0.02) X 10-^ (115) 

Aahad(T(25-65)) = (0.62 ± 0.03) X 10"^ (116) 



3.11 Inclusive hadronic data contribution {^/s < 11.09 GeV) 

We use four data sets below 2 GeV jHlj-IEZI, and twelve data sets above 2 GeV |HH]^|nH] (see 
Fig.lTBj). Below 2 GeV, we correct for the unaccounted modes. Namely, we add the contributions 
from the u!{^ 7r°7)7r'^ and Kg{^ 2tt^)K^tt^ channels to the experimentally observed i?-ratio, 
since the final states of these channels consist only of electrically neutral particles, which are 
hard to see experimentally. They shift the R values by roughly 1%, depending on ^/s. In 
addition we correct some experiments for the contributions from missing two-body final states, 
as discussed at the end of Section 2. We have also checked that corrections for 7 — Z interference 
effects are completely negligible in the energy range below 11.09 GeV where we use data. 
The contributions to the muon g ■- 2 and Atthad are, from 1.43 < y/s < 2 GeV, 

a^(inclusive, ^/s < 2 GeV) = (31.91 ± 2.42) x 10"^°, (117) 
Aohad (inclusive, < 2 GeV) = (10.78 ± 0.81) x 10-^ (118) 

and from 2 < < 11.09 GeV, 

a^(inclusive, 2 < < 11.09 GeV) = (42.05 ± 1.14) x 10"^°, (119) 
Aahad (inclusive, 2 < < 11.09 GeV) = (81.97 ± 1.53) x 10'^ (120) 
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respectively. 



3.12 Inclusive pQCD contribution (^s > 11.09 GeV) 

Above 11 GeV we use perturbative QCD to evaluate the contributions to a^'^^ and ^a{M'^). 
We incorporate 0{a.^g) massless quark contributions, and the massive quark contribu- 

tions |107| I124L I125[ I126[ I127j . We have checked that our code agrees very well with the code 
rhad written by Harlander and Steinhauser |128j . As input parameters, we use 

a5(M|) = 0.1172 ±0.002, = 174.3 ± 5.1 GeV, m;, = 4.85 ± 0.25 GeV, (121) 

and allow for an uncertainty in the renormalization scale of ^/i/2 < n < 2^/s. Here rrit and rrib 
are the pole masses of the top and bottom quark. We obtain 

a^(pQCD, > 11.09 GeV) = (2.11 ± 0.00) x 10~^°, (122) 

where the uncertainty from as{M'^) is dominant, which is less than 1 x 10^^^. Similarly, for 
Acthad we find 

Aahad(pQCD,v^> 11.09 GeV) = (125.32 ± 0.14 ± 0.02 ± 0.01) x 10^^ (123) 

= (125.32 ± 0.15) X 10"^ (124) 

where the first error comes from the uncertainty in as{M'^), the second from the renormalization 
scale n, and the third from that on the mass of the bottom quark. 

3.13 Total contribution to the dispersion integrals 

To summarize. Table shows the values obtained for a^^^'^'LO g^j^^ Acthad? as well as showing the 
contributions of the individual channels. Summing all the contributions we obtain 

aj^<^'^°(incl.) = (692.38 ± 5.88cxp) X 10"^°, (125) 
^had,LO(^g^pl ^) = (696.15 ± 5.68exp) X 10^1°, (126) 

where "inch" means that we have used the inclusive data sets for 1.43 < y/s < 2 GeV, while 
"excl." means that we used the exclusive data at the same interval, "exp." means that the 
errors are from the experimental uncertainty. The corresponding results for Aa^ad are 

Aahad(incl.) = (275.52 ± 1.85exp) x 10-^ (127) 
Aahad(excl.) = (276.90 ± 1.77cxp) x 10"^ (128) 

We see that using the sum of the data for exclusive channels to determine R{s), in the in- 
termediate energy interval 1.43 < y/s < 2 GeV, yields values for a^^"^'^^ and Acthad which 
significantly exceed the values obtained using the inclusive data for R{s). The mean values 
differ by about 2/3 of the total experimental error. In Fig. El we show the hadronic R ratio as 
a function of ^/s. A careful inspection of the figure shows the discrepancy between the inclusive 
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and exclusive data sets in the interval 1.43 < i/i < 2 GeV, see Fig. |3] The contribution from 
this region alone is 



ap^°(1.43 < < 2 GeV, inch) = (31.91 ± 2.42e,p) x 10"^°, (129) 
^had,LO^^_43 < ^ < 2 GeV, excl.) = (35.68 ± l.Tlcxp) x 10~^°, (130) 

and for Aa^ad, 

Aahad(l-43 < < 2 GeV, inch) = (10.78 ± 0.81exp) x 10-^ (131) 
A«had(l-43 < < 2 GeV, excl.) = (12.17 ± 0.59cxp) x 10"^ (132) 

In the next Section we introduce QCD sum rules that are able to determine which choice of 
R{s) is consistent. We find that the sum rules strongly favour the use of the inclusive data in 
the above intermediate energy interval. 

Table U\ shows the breakdown of the contributions versus energy. It is also useful to show 
the breakdown visually in terms of 'pie' diagrams. The 'pie' diagrams on the left-hand 
side of Fig. show the fraction of the total contributions to aJ^^'d-LO Aahad coming from 
various energy intervals of the dispersion integrals (jlj) and The plots on the right-hand-side 
indicate the fractional contributions to the square of the total error, including the error due to 
the treatment of radiative corrections. The values shown for a^^'^'^'^ in these plots correspond 
to using the inclusive data in the intermediate energy interval. 

In Section 13 we use the value of a^^'^'^'^, along with the QED, weak and other hadronic 
contributions, to predict the value of (7 — 2 of the muon. In Section |H1 we use the value of 
Aahad(^I) to predict the value of the QED coupling on the Z pole, a(M|). 



4 Resolution of the ambiguity: QCD sum rules 

To decide between the exclusive and inclusive data in the energy range 1.43 < y/s < 2 GeV 
(see Fig. Ej), we make use of QCD sum rules |129j . see also the review |130j . The sum rules are 
based on the analyticity of the vacuum polarization function n(g^), from which it follows that 
a relation of the form 



/ dsR{s)f{s) = / dsD{s)g{s) (133) 



must be satisfied for a non-singular function f{s). C is a circular contour of radius sq and g{s) 
is a known function once f{s) is given. The lower limit of integration, Sth, is 4m^, except for a 
small e^e~ — >■ 7r''7 contribution. D{s) is the Adler D function, 

Dis) = -UTT^s^i^^] , where R(s) = — Imn(s). (134) 
as \ s J s 

Provided that sq is chosen sufficiently large for D{s) to be evaluated from QCD, the sum rules 
allow consistency checks of the behaviour of the data for R{s) for s < sq. Indeed, by choosing 
an appropriate form of the function f[s) we can highlight the average behaviour of R{s) over 
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Figure 17: The hadronic R ratio as a function of ^/s. Note that the values of R obtained from 
the sum of the exclusive channels and from the inclusive data overlap in the region 1.4 ^ ^ 2 
GeV. 
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value (error)^ 




Figure 18: The pie diagrams in the left- and right-hand columns show the fractions of the 
total contributions and (errors)^, respectively, coming from various energy intervals in the 
dispersion integrals ^ and (0). The pie diagrams for the LO hadronic contribution to 
g — 2, shown in the first row, correspond to sub-contributions with energy boundaries at 
mjr, 0.6, 0.9, 1.4, 2 GeV and oo, whereas for the hadronic contribution to the QED coupling, 
shown in the second row, the boundaries are at m^r, 0.6, 0.9, 1.4, 2, 4, 11.09 GeV and oo. In the 
(error) ^ pie diagrams we also included the (error) ^ arising from the treatment of the radiative 
corrections to the data. 
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a particular energy domain. To be specific, we take Sq just below the open charm threshold 
(say y/sQ = 3.7 GeV) and choose forms for f{s) which emphasize the most ambiguous range 
(1.5 ^ ^/s ^ 2 GeV) of R{s), so that the discriminating power of the sum rules is maximized. 
We therefore use the three flavour (n/ = 3) QCD expressions for D{s), and omit the J/ip and 
ipi^S) cc resonance contributions to R{s). 

To evaluate the function D{s) from QCD, it is convenient to express it as the sum of three 
contributions, 

Dis) = Do{s) + D^{s) + Dr,p{s), (135) 

where Dq is the 0(a|) massless, three- flavour QCD prediction, is the (small) quark mass 
correction and D^p is a (very small) contribution estimated using knowledge of the condensates. 
Dq is given by |124j 

D„(-.) = aEQj {l + ^ + 4 (^)' + 4 (^)' + O (a|(.))| . (136) 



with 



di = 1.9857 - 0.1153n/, (137) 
d, = d2 + ^ fwith/5o = ll-^V (138) 
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d2 = -6.6368 - 1.2001nj - 0.0052nJ - 1.2395 '^-^ ^''^ , (139) 

where the sum / runs over u, d and s flavours. Q/ is the electric charge of quark /, which takes 
the values 2/3, —1/3, and —1/3 for u,d and s, respectively. The quark mass correction Dm 
reads |131j 

D^{-s) = -3 E fe + 28^ + (294.8 - 12.3n,) ( ^] . (140) 



/ 



We take the MS s-quark mass at 2 GeV m^(4 GeV^) to be 120 ± 40 MeV, and we neglect the 
u and d quark masses. The contribution from condensates, -Dnp, is given by 



+ 87r2 (l - ^s{s) \ {mfqfqf) 327r^ asjs) ^ {nikg^qk) 
\ TT J 27 TT ^ s2 

+ 12^^i^ + 16n'^ \ , (141) 



where, following |132j . we take 

'as 



GG) = 0.037 ±0.019 (GeV 

2_2 



{m-ss) = -f^mj,. (142) 
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Here /^r — 92 MeV is the pion decay constant, and mx is the kaon mass. As we will see 
later, the quark mass corrections and the condensate contributions are very tiny — typically at 
most a few percent of the whole QCD contribution. Hence we neglect the higher dimensional 
condensates, {Oq) and (Cs)- 

As for the weight function /(s), we take it to be of the form (1 — s/so)™'(s/so)" with 
n + m = 0, 1 or 2. For these six choices of /(s), the function g{s) may be readily evaluated, 
and the sum rules, ()133p . become 




We evaluate each of these sum rules for ^/sq = 3.7 GeV using the clustered data values 
of R{s) of Section 2 on the l.h.s., and QCD for D{s) (with ag = 0.1172 ± 0.0020 [Elj) on 
the r.h.s. We find, as anticipated, that the sum rules with m = and n = 1 or 2 have very 
small contributions from the disputed 1.43 — 2 GeV region. Indeed, this region contributes only 
about 5% and 2%, respectively, of the total contribution to the l.h.s. of ()144|1 and ()146|1 . They 
emphasize the region s ^ sq and so essentially test data against perturbative QCD in this small 
domain. They are not useful for our purpose. The results for the remaining four sum rules are 
shown by the numbers in brackets in Fig. ^| For this choice of sq, the sum rules with m = 1 
or 2 and n = are found to maximize the fractional contribution to the sum rule coming from 
the 1.43 < < 2 GeV interval. These two sum rules clearly favour the inclusive over the 
exclusive data. 

The comparison between the data and QCD can be translated into another form. We can 
treat a5(M|) as a free parameter, and calculate the value which makes the r.h.s. of a sum rule 
exactly balance the l.h.s. The results are shown in Fig. El We can see that in this comparison 
the determination from the inclusive data is more consistent with the world-average value, 
as{M^) = 0.1172 ± 0.0020 [Tllij . 

For illustration, we show in Table |H1 a detailed breakdown of the contributions to both 
sides of the sum rule for the cases of m = 2, n = and m = n = 0. If we compare the 
breakdown of the contribution from the data in both cases, we can see that the weight function 
f{s) = (1 — s/sq)"^ highlights the most ambiguous region of R{s) very well. When we look into 
the breakdown in the QCD part, we can see that the QCD contribution is dominated by the 
massless part. 

We repeated the sum rule analysis for ^/sq = 3.0 GeV, see Fig. UHl The lower value of 
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So means that more weight is given to the disputed 1.43 — 2 GeV region. Taken together, we 
see that the sum rules strongly favour the behaviour of R{s) from the inclusive measurements. 
Indeed, the overall consistency in this case is remarkable. This result can also be clearly seen 
from Fig. ^[ which compares the world average value of as{M^) with the predictions of the 
individual sum rules for, first ^/sQ = 3.7 GeV and then for ^/sq = 3 GeV. Again the consistency 
with the inclusive measurements of R{s) is apparent. 



1 I I I I I I I rr| I I I I I I I I I I I I I I I I 
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Figure 19: The QCD sum rule predictions for as{M'^) compared with the world average 
value |lU4j . The results for the four sum rules for two values of ^/sq are shown. In each 
case we show results for the inclusive and the exclusive measurement of R{s) in the intermedi- 
ate energy region. We also give in brackets the fractional contribution to the sum rule coming 
from the 1.43 < ^/s < 2 GeV interval. 

The same conclusion with regard to the resolution of the inclusive/exclusive ambiguity in 
the 1.43 < ^/s < 2 GeV interval was reached in an independent analysis |133j . 

In an attempt to understand the origin of the discrepancy, we have studied the effect of 
possibly missing (purely neutral) modes in the inclusive data, but found that these cannot 
explain the difference. One should, however, keep in mind that the precision of both the (old) 
inclusive and the exclusive data in this energy regime is quite poor. We expect that future 
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measurements at 5-factories (via radiative return) and at the upgraded machine VEPP-2000 
in Novosibirsk will improve the situation in the future. 

5 Comparison with other predictions oi g — 2 

Fig. 1201 shows other determinations of a^^'^'^'^, together with the values (HMNT(03)) obtained 
in this work. The values listed below the first dashed line incorporate the new more precise 
data on e^e~ — > tt^tt" into the analysis. These data play a dominant role, and, as can 
be seen from the Figure, significantly decrease the value of a^^'^. However, very recently, the 
CMD-2 Collaboration have re-analysed their data and found that they should be increased by 
approximately 2%, depending on ^/s. The new data [TU] are included in our analysis. Inspection 
of Fig. EO shows that the re-analysis of the CMD-2 data has led to an increase of a^^'^'^'^ x 10^° 
by about 10. 
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Figure 20: Recent evaluations of aj'^'^'^o H 1121 lEEl CHSl CESl IMl M^- The entries 
below the first dashed line include the new CMD-2 tt+tt" data JT], and the values below the 
second dashed line include the re-analysed CMD-2 tt+tt^ data PH] . 



The entries denoted by "DEHZ (r)" also used information from hadronic r decays [3 El, 
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which through CVC give independent information on the e~^e~ 27r and An channels for 
^/s < 171^-. The apparent discrepancy between the prediction from this analysis and the pure 
e~^e~ analyses is not yet totally understood, however see the remarks in the Introduction. 

5.1 Comparison with the DEHZ evaluation 

It is particularly informative to compare the individual contributions to a^'^'^'^'~' obtained in the 
present analysis with those listed in the recent study of Davier et al. (DEIIZ03) [Sj, which used 
essentially the same e~^e~ — * hadrons data. Such a comparison highlights regions of uncertainty, 
and indicates areas where further data and study could significantly improve the theoretical 
determination of g — 2. DEHZ provided a detailed breakdown of their contributions to a^^'^'^'^, 
and so, to facilitate the comparison, we have broken down our contributions into the energy 
intervals that they use. Table El shows the two sets of contributions of the individual e~^e~ 
channels to dispersion relation in the crucial low energy region with ^/s < 1.8 GeV. 

The last column of Table El shows the discrepancy between the two analyses. The biggest 
difference occurs in the tt+tt^ channel, which gives the main contribution to a^^'^'^'^, and the 
improvement in the SM prediction essentially comes from the recent higher precision CMD- 
2 data in the region 0.6 < ^/s < 0.9 GeV (see the remarks in Chapter 13. 2j) . We find that 
this difference, 2.6 x 10^^°, appears to come from the region just above the vr+vr" threshold, 
especially in the region y/s ~ 0.4 GeV, see Fig. |^ The figure shows the vt+tt^ contribution 
plotted in such a way that the area under the curves (or data band) gives the contribution to 
dispersion relation for a^^'^'^'^. To determine the low energy vr+vr" contribution, DEHZ |2| 
first perform a three-parameter fit to vr^vr" data for ^/s < 0.6 GeV, and obtain the dashed 
curve in Fig. 1^ This is then used to compute the tt+tt" contribution of (58.04 ±2.06) x 10~^° 
for y/s < 0.5 GeV. They do not use either the NA7^^ [201 or the preliminary CMD-2 data. On 
the other hand we use the chiral description |136j . shown by the continuous curve, only as far 
as y/s = 0.32 GeV; and then use the band obtained from our clustered data, which include 
data from OLYA TOF NA7 |2D|, CMD jHl and DM1 in this energy region. In 
this way we obtain a vr+vr" contribution for ^/s < 0.5 GeV of (55.7 ± 1.9) x 10~^°. We also 
show on Fig. |^the preliminary CMD-2 data, obtained from Fig. 3 of Ref. |134j . These data 
were used in neither analysis, but do seem to favour the lower tt+tt^ contribution. It is also 
interesting to note that DEHZ |21IS| obtain the low value of (56.0 ± 1.6) x 10~^° if r decay and 
CVC are used in this region. 

Other significant, with respect to the errors, discrepancies arise in the tt^j + 777 and the 
KgKl channels, where the treatment is different: DEHZ integrate over Breit-Wigner resonance 
parametrizations (assuming that the KK channels are saturated by the decay), while we are 
integrating the available data in these channels directly. In our method there is no danger to 
omit or double-count interference effects and resonance contributions from tails still present at 
continuum energies, and the error estimate is straight forward. As a check, we made fits to 
Breit-Wigner-type resonance forms and studied the possibility that trapezoidal integration of 
concave structures overestimate the resonance contributions. We found the possible effects are 

^^However, recently it has turned out that earlier worries about a systematic bias in the NAT data as mentioned 
in [2] are not justified and that there is no reason to neglect these important data |135| . 
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Figure 21: The tt+tt^ data just above threshold, plotted so that the area gives the contribution 
to dispersion relation (ji^ for a^^'^'^'^. The dashed curve is used by DEHZ whereas the 

continuous curve up to y/s = 0.32 GeV (s = 0.10 GeV^) and data band are used in this analysis; 
see text. We also show, but do not use, the preliminary low energy CMD-2 data, which were 
read off Fig. 3 of |134j . These points, particularly the first, are subject to 'reading-off ' errors. 
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negligible compared to the uncertainties in the parametrization coming from poor quality data 
in the tail regions. The one exception is the — >■ KgK^ contribution. Here the lack of data in 
certain regions of the resonance tails (see Fig. IT3|) has caused us to increase the uncertainty on 
this contribution to a^, see Section ESI 

Apart from these channels, it is only the two four-pion channels which show uncomfortably 
large and relevant discrepancies. Here, the data input is different between DEHZ and our 
analysis. We use, in addition to DEHZ, also data from 772 ISHIESI and ORSAY-ACO jSSl for 
both An channels, and data from MSN jSU] and two more data sets from CMD-2 [UZl EH] for 
the 71^71^71^71^ channel. However, it should be noted, that the available data are not entirely 
consistent, a fact reflected in the poor x^j^/d.o.f. of our fits resulting in the need of error 
infiation^^. Clearly, in these channels, new and better data is required. As mentioned already 
in Section Em the situation is expected to improve as soon as the announced re-analysis from 
CMD-2 will become available. 

There are no data available for some of the exclusive channels. Their contribution to the 
dispersion relation is computed using isospin relations. The corresponding entries in Table El 
have been marked by the word "isospin" . 



5.2 Possible contribution of the (j(600) resonance to g — 2 

This subsection is motivated by the claim |137j that the isosinglet scalar boson^° (t(600) can 
have a non-negligible contribution to the muon g — 2. Here, we evaluate its contribution and 
find that it is at most of order 0.1 x 10~^°. This is negligible as compared to the uncertainty of 
the hadronic vacuum polarization contribution of 6 x 10~^°, and hence we can safely neglect it. 

The argument presented in Ref. |137j is twofold. First a may contribute to the muon g — 2 
through unaccounted decay modes of the narrow spin 1 resonances into the channel. The 
second possibility, considered in |137j . is that a may contribute directly to the muon g — 2 
through its coupling to the muon pair. We estimate the two contributions below. 

In the zero-width limit, narrow spin 1 resonances, V, contribute to the muon g — 2 as 

al = i3/7i)K{m^y)T{V ^ ee)/mv, (149) 
where K{my) is the kernel function (j45j) at s = my. We find, for example^^, 

a'^^ = 391 X 10"^°, (150) 

aj = 39 X 10"^°, (151) 

^0(1.68) _ 3 4 ^ ^^52) 

where, in ()152|) . we have used r(0(1.68) — > ee) = 0.48 keV j.104^ to give a rough estimate. 
If the decays V —>■ a'-f of the above vector bosons escape detection, a fraction of the above 

-'^^If for a given channel Xmin/d-o.f. > 1.2, then we enlarge the error by y^Xmin/'^-'^-^-- This was necessary for 
three channels, see Section ITU 

2°(t(600) is denoted by /o(600) in the Review of Particle Physics [TUl) . 

^^We take vector mesons, V, which, according to |137| . may have significant contributions. 
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contributions np to B{V ^ a^) may have been missed. On the other hand we find that 99.8% 
of decays has been accounted for in the five decay channels exphcitly included in our analysis 
hence B{(f) — * (77) < 0.002. This severely constrains the (T77 coupling. Hence we can use the 
Vector Meson Dominance (VMD) approximation to show that the other branching fractions 
satisfy B{u ^7) < 7.2 x 10~5 and 5(0(1.68) ^ ^7) < 3.5 x 10-^ see Appendix B. By 
inserting these constraints into the estimates (jl50j) - (jl52p . we find that the unaccounted ^ (T7 
contributions to (7 — 2 of the muon are less than (2.8, 7.8, 0.01) x 10~^^ for V = oj, 0, 0(1.68) 
respectively, assuming = 600 MeV. These estimates are much smaller than those presented 
in |137j . It is clear that the total contribution of unaccounted (T7 modes through narrow 
resonance decays is negligibly small. It is also worth pointing out here that the unaccounted 
fraction 0.2% of the contribution (7.8 x 10^^^ above) has been taken into account in our 
analysis, whether it is — > (17 or not. 

We now turn to the a contribution to the muon g — 2 through its direct coupling to a muon- 
pair. To evaluate this, it is essential to estimate the magnitude of the a/i/i coupling. Since the 
coupling through the a-Higgs boson mixing is negligibly small, the leading contribution should 
come from two-photon exchange. In this regard, the effective coupling strength should be of 
the same order as that of the rj isoscalar pseudoscalar meson, which should also be dominated 
by the two-photon exchange. By using the observed width V{rj /i/i), and by neglecting the 
form factor suppression, we find that the point-like r] contribution to the muon g — 2 is 

a;^ = -3 X 10"^^ (153) 

which is negligibly small. It follows that this implies that is also negligibly small, see ()164p 
below. However, the discussion can be made far more general. It is presented in the next 
Section. 

6 Internal light-by-light contributions 

In this section we present a very primitive discussion of the hadronic contribution to the internal 
light-by-light amplitudes, motivated by the study of the direct a and i] contributions to the 
muon g — 2. 

The meaning of 'internal' can be seen from Fig. E21 We call the diagram on the right 
'internal' to distinguish it from the left diagram, that is the familiar light-by- light contribution 
which, here, we call 'external'. We should note that the external light-by-light diagram is of 
0{a^) and the internal light-by-light diagram is an 0(q;^) contribution. 

6.1 Internal meson contributions 

Just like the external light-by-light amplitude is dominated by a single pseudoscalar meson 
contribution |138j - [TiT] . it is likely that the hadronic contribution to the internal light-by-light 
amplitudes is dominated by a single meson exchange contribution. 

In general, we can estimate the internal contribution to from arbitrary scalar and pseu- 
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external internal 

Figure 22: External and internal light-by-light contributions to g — 2. The former is an 0{a^) 
and the latter is an 0{a^) contribution. In this paper we compute the contribution of the 
internal diagram. In Section 16.11 we take the shaded blob to be scalar (a) or pseudoscalar 
(tt",?]) mesons, whereas in Section we take it to be a light {u,d,s) quark loop, using the 
result for lepton (e, /i) loops clS db guide. 



doscalar mesons. Using the effective coupling 



we find [132 
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where r = m^/m|, with h = S and P in ()155p and ()156|) respectively. The scalar contribution 
is positive definite and the pseudoscalar contribution is negative definite. In the large mass 
limit (r ^ 1) we have 
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Further, in the parity-doublet limit oi gs = gp and nis = nip, the leading terms cancel 
and only a tiny positive contribution remains. The effective couplings in ()154|) can be extracted 
from the leptonic widths 
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where n = 3, 1 for h = S,P respectively. 

Let us estimate the pseudoscalar 7r° contribution. We use 



e e 



5 X 10"^ eV. 



(160) 



After we allow for the helicity suppression factor of nie/m^ for the n^ee coupling, this gives a 
Ti^fifi coupling 
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and hence, from ()156|) . a contribution 

af ~ -6 X 10-^\ (162) 

Although this contribution is not completely negligible, we expect a form factor suppression 
of the effective couplings and so the pion structure effects should suppress the magnitude 
significantly. 

In the scalar sector, we do not find a particle with significant leptonic width. Although the 
a leptonic width is unknown, we find no reason to expect that its coupling is bigger than the 
r]lili coupling. If we use 

^ ^ i = ^(^^ ^ ^^^'^ ^ 1.4 X 10-- (163) 
Svr Svr m^Wl - Amf^/m^ 

we find that r((T ^ /^^/^ ) — 7 x lO^'^eV, and hence 

a)^ = 7 X 10-^^ (164) 

for m„ = 600 MeV. Again we should expect form-factor suppressions. Because pseudoscalar 
mesons are lighter than the scalars, there is a tendency that the total contribution is negative 
rather than positive. 



6.2 Internal lepton or quark contributions 

The internal light-by-light scattering contributions in the 4-loop order have been evaluated in 
QED. The electron-loop contribution is 



a 



,int. 1-b-l/ 



e-loop) = -4.43243(58) f-V ~ 1.29 x 10"^°, (165) 



TT 



whereas the muon-loop contribution is 

~ -0.29 X 10-'°. (166) 
The /i-loop contributions to the electron anomalous moment has also been estimated |145j 

■ ^-^"V-loop) = -0.000184(14) (^^^ . (167) 

If we interpolate between p66|) and p67|) by assuming the form {m'j^/mf)[A\n{m'f /m'j^) + B], 
we obtain the estimate 



a}^'- ^-^-'(/-loop) ^ - 



m 



0.65 In ^ + 1 







mf 





'168) 



which may be valid for an arbitrary lepton mass in the range < rrii < m?^/me ~ 20 GeV. 
For a r-loop internal light-by-light contribution to a^, the relation (|168|) gives 



int. i-b-i^^_iQop) ^ _o.0165 (^j , (169) 
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which agrees with the actual numerical result 



int. 1-b-l, 
''/I ' 



r-loop) = -0.01570(49) - 



a 



TT 



(170) 



within 10%. We can now estimate the hadronic contribution by using the constituent quark 
model 

2 



int. 1-b-l/ 



U, d, S-Ioop) — —77 

O 



r 2 1 

TTl 


^ / 


0.65 In ^ + 1 




m2 





(171) 



where we use rriu = rud = mg = rrin to set the scale, and where | = 3 ( (|) + 2 (\") ) is the 



charge factor. For rriq = 300 MeV, ()17Hl gives 

^'^'Ku, d, s-loop) ~ -6 X 10-i2_ 



(172) 



6.3 Quark loop estimates of the hadronic Ught-by-hght contribu- 
tions 

If the same massive quark loop estimate is made for the 3-loop (external) hadronic light-by-light 
scattering contribution, it is found that |146j 

a;f*- '-^-1(m, d, s-loop) ~ - X 0.615 ('^] (-Y -6x 10"^°. (173) 

As we shall see later, this estimate is in reasonably good agreement with the present estimate 
of the total contribution of (8 ± 4) x 10"^° of ()192|) . and of its sign. 

The above well-known result has been regarded as an accident, because in the small quark 
mass limit the quark-loop contribution to the external light-by-light amplitude diverges. The 
light-meson contributions could only be estimated by adopting the effective light-meson de- 
scription of low-energy QCD. Although the same may well apply for the internal light-by-light 
amplitudes, we note here that the quark-loop contributions to the internal light-by- light ampli- 
tudes remain finite in the massless quark limit because of the cancellation of mass singularities 
|147t I148j . We find no strong reason to discredit the order of magnitude estimate based on 
f)172|) against the successful one of ()173p for the external light-by-light amplitudes. Although 
the point-like vr contribution of ()162|) is a factor of ten larger than the estimate ()172|) . the cor- 
responding point-like vr contribution to the external light-by-light amplitudes diverges. We can 
expect that the form factor suppression of the effective vertices should significantly reduce its 
contribution. Also, since these mesons are lighter than the scalar mesons, we expect the sign of 
the total meson contribution to be negative, in agreement with the quark loop estimate of p72|) . 
In conclusion, we use fll72j) to estimate that the hadronic internal light-by-light contribution is 
given by 

■ ^-^-^(hadrons) = -(0.6 ± 0.6) x 10"", (174) 
which is totally negligible. We do not take this contribution into account in our final results. 



52 



7 Calculation of ajj^^ and g — 2 oi the muon 
7.1 Results on a^^^'^^ 

We calculated the LO hadronic contribution aJ^a,d,LO Sectional We found 

^had,LO _ (692.4 ± 5.9exp ± 1.4rad,vp ± I-Wsr) X 10-^° (175) 
= (692.4 ± 5.9exp ± 2.4rad) X 10-1°, (176) 

where the first error comes from the systematic and statistic errors in the hadronic data which 
we included in the clustering algorithm, and the second error is from the uncertainties in the 
radiative correction in the experimental data. Below we explain this in more detail. 

We add the VP error from the experiments and the narrow resonances linearly. Out of 
1.4 X IQ-i", 1.2 X 10"^° is from the data, and 0.2 x 10~^° is from the narrow resonances. For 
the errors from the final state radiation we assign 1.9 x 10^^*^, which is the sum of the errors, 

^^fsr.TT+TT- ^ Q gg X ^Q-IO^ ^^lsr,K+K- ^ g_42 ^ 10~W ^-^^ ^^fer, other excl. ^ Q_g^ ^ 

added the errors from the VP and the FSR in quadrature, which is the second error in ()176j) . 



7.2 Calculation of the NLO hadronic contributions to — 2 of the 
muon 

In this subsection we update the computation of the NLO hadronic contribution to (7 — 2 of the 
muon. It proceeds in a similar way to that for the LO contribution, but now the kernel of the 
dispersion relation is a little more complicated. There are three types of NLO contributions, 
which were denoted (2a), (2b) and (2c) by Krause |149j : (2a) consists of the diagrams which 
contain one hadronic bubble and which do not involve leptons other than the muon, (2b) is 
the diagram which has one hadronic bubble and one electron (or tau) loop, and, finally, (2c) is 
the diagram which has two hadronic bubbles. The three different classes of NLO contributions 
correspond to the diagrams which are denoted (a,b,c) respectively in Fig. 1221 
The contributions from (2a), (2b), and (2c) can be written as 

/ dsaUs)K^"^\s), (177) 



had,NLO(2a) 



47r4 
a 

4^ 



Sth 

00 



a™^^^^ = r'dsaUs)K^''\s), (178) 

'Sth 



1 poo roo 

«r''''°^'^^ = TfiV/ / ds'aUs)aUs')K^"'\s,s'), (179) 
167r^a J Jsth 

where the analytic expressions for if (2") , if (2^) and K^'^^^ are given in Ref . |149j . We use the 
clustered data for the cross section for e^e" hadrons, a^^^^ of (j21), with the choice of inclusive 
data in the regime above 1.43 GeV to compute the contributions of the three different classes 
of NLO diagrams. We find 

^had,NLO(2a) _ (-20. 73 ± 0. 18cxp ± 0.07rad) X 10-^°, (180) 

^had,NLO(2b) _ (^Q gQ ^ Q Qg^^^ ^ 0.04rad) X 10"^°, (181) 
^had,NLO(2c) _ (0.34 ± O.Olexp ± O.OOrad) X 10-1°, (182) 
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had. had 





(b) i (c) i 

f ^^^had. had^^had. 



Figure 23: The three classes of diagrams (a,b,c) which contribute to a^ad.NLO^ Qjg^gg (^^^^ contains 
the first five diagrams. In the class (b) diagram, / = e or r, but not fi. Mirror counterparts 
and diagrams with an interchange between the massless photon and the "massive photon" 
propagators should be understood. 



where we have assigned the uncertainty from the radiative correction similarly to the LO 
hadronic contribution. When summed, 



^had,NLO _ (-9.79±0.09exp±0.03,ad) X 10^^ 

which may be compared to the original calculation of Krause |149j . 

^had,NLO ^ (_21. 1(0.5) + 10.7(0.2) + 0.27(0.01)) x 10"^° = -10.1(0.6) x 10"^°. 



;i83) 



In ()183p we added the error linearly with an opposite relative sign since the errors in (2a) 
and (2b) are nearly 100% correlated in the opposite directions. Hence the total error is the 
difference of the two. In combining the errors we neglected the errors from (2c) since it is 
negligibly small compared to the other errors. 

Note that the contribution of diagram (2c) does not agree with the result given by Krause, 
when account is taken of the small error on this contribution. We have therefore performed two 
checks of our numerical programme. First we replaced the two hadronic blobs of the diagram 
(2c) with two muon loops, since the contribution from such a diagram is known analytically 
|150j as a part of the QED contribution. It is 



a^(two muon loops along one photon propagator) 
943 

.tt) V 324 ~ 
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^C(2) + 3C(3) 



184) 



= (^^) X 0.002558... (185) 
= 0.3206... X 10"^°. (186) 

Our programme reproduced 0.321 x 10^^°, which agrees with p86|) within an accuracy of 10~^^, 
which is the accuracy of the calculation throughout this paper. 

As a second check, we have taken R{s) to be a step function. In the first line of Eq. (13) 
of the paper by Krause, the contribution from the diagram (2c) is written as a triple integral 
over s, s' and x, where s and s' are "mass- squared" of the hadronic blobs, and x is a Feynman 
parameter. By explicitly integrating over x, Krause obtained the second line of Eq. (13), which 
is a double integral over s and s'. We are using this expression to integrate over the hadronic 
data. If R{s) is a constant, we can explicitly integrate over s and s', instead of x. Then we 
are left with one dimensional integral over x, which is much more tractable than the double 
integral over s and s'. We compared the result obtained from this integral over x with the 
double integral over s and s'. Below are the numerical results. 

When R{s) is a constant (more rigorously, when R{s) is a step function with R{s) = 1 for 
s > 4m^, otherwise R{s) = 0), the result from the double integral is 

= 0.21 X 10^^°, (187) 

(which has only two significant digits) and the result from the integral over x is 

= 0.2109... X 10"^°. (188) 

The agreement is very good. From the above two checks we believe our result for diagram (2c) 
is correct. 



7.3 Hadronic contribution to (/ — 2 of the muon 

The hadronic contribution a^^'^ has been divided into three pieces, 

^had ^had,LO i ^had,NLO i ^had,l-b-l /'1Qn^ 

= ' ' + ' • (189) 

The lowest-order (vacuum polarisation) hadronic contribution, a^^'^'^'^, was calculated in Sec- 
tional There we found 



a 



had,LO ^ ^gg2,4 ± 5_9^^p ± 2.4rad) X 10"^°, (190) 



where we have used the QCD sum rule analysis to resolve the discrepancy in favour of the 
inclusive e~^e~ — > hadrons data in the region 1.4 < y/s < 2 GeV. The value of the next-to- 
leading order hadronic contribution, a^^'^'^^'~', was updated by the calculation described in the 
previous subsection. We obtained 

^had,NLO ^ (-9.79 ± 0.09exp ± 0.03,ad) X 10-1°. (191) 
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Finally, we must include the hadronic light- by-light scattering contribution a^^'^'^'^'^. It has 
attracted much study. Recent re-evaluations can be found, for example, in Refs. |151j - P^Bj . 
Here we take the representative value^^ 

as given in Ref. |157j . From Eqs. ()176|) . (jl91|) and ()192p . we can see that a^^'^'^^ has the largest 
uncertainty, although the uncertainty in the light-by-light contribution also large. 

When we combine all the three contributions to the hadronic contribution, we find 

aj'^d = (690.6 ± 7.4) x 10"^°. (193) 

To calculate the number above, we first added the uncertainties associated with the LO and 
NLO diagrams linearly, and then added the uncertainty in the light-by-light contribution 
quadratically. We did so since the errors in the LO and the NLO contributions are nearly 
100% correlated. 



7.4 SM prediction of — 2 of the muon 

The SM value of the anomalous magnetic moment of the muon, a^, may be written as the sum 
of three terms. 

The QED contribution, a^^^, has been calculated up to and including estimates of the 5-loop 
contribution, see reviews |144( 11581 1159| 1160] , 

aj^° = 116 584 703.5(2.8) x 10"". (195) 

This value |160j includes the recent update from |144j . In comparison with the experimental 
error in Eq. (^J), and the error of the hadronic contribution, the uncertainty in a^^^ is much 
less important than the other sources of uncertainty. The electroweak contribution aj^^ is 
calculated through second order to be |161j - P^ij 

qEW ^ ^54(2) X 10-11. (196) 

Here we quote the result of |164j . Although some discrepancies on conceptual questions remain, 
this result agrees numerically with the one of |163j , and here again the error is negligibly small. 

Summing up the SM contributions to a^, as given in ()193|) . ()195|) and ()196|) . we conclude 
that 

af^ = (11659176.3 ± 7.4) x 10"i°, (197) 

which is 26.7 x 10"^'^ (2-4cr) below the world average experimental measurement. If, on the 
other hand, we were to take, instead of (jl90|) . the value of ah^^.LO obtained using the sum of 
the exclusive data in the interval 1.43 < ^/s < 2 GeV, then we would find = (11659180.1 ± 
7.4) X 10^1°, which is 22.9 x 10"^° (2. la) below a^^'P. The above values of a^^^ is compared 
with other determinations in Fig. |211 
^^However, see the note added in proof. 
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Figure 24: Recent evaluations of and the current world- average of the measured value 
(shown as a band). The band corresponds to a l-cr range. The final values, HMNT(03), are 
the predictions of this work, and include the recently re-analysed CMD-2 tt+tt" data in our 
analyses. 

8 Determination of aQED(^|) 

As mentioned in the Introduction, the value of the QED coupling at the Z boson mass is 
the least well known of the three input parameters (G^, Mz and a{M'^)) which are the three 
most fundamental inputs of the standard electroweak model. Its uncertainty is therefore the 
limiting factor for precision electroweak physics. It is clearly important to determine a(M|) as 
accurately as possible. 

The value of a(M|) is obtained from jlU4j 

a-^ = a{<d)-^ = 137.03599976(50) (198) 

using the relation 

aisy' = (l - Aaiep(s) - Aa^t^s) - Aa*°P(s)) a-\ (199) 
where the leptonic contribution to the running of a is known to three loops |169j . 

Aaiep(M|) = 0.03149769. (200) 
The evaluation of the hadronic contribution, Aa^^^^{M'^), is described below. 
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8.1 The hadronic contribution to the running of a up to s = M| 

It is conventional to determine the contribution from 5 quark flavours, Aa^^g^, and to include 
the contribution of the sixth flavour |17Uj , 

Aa*°P(M|) = -0.000070(05), (201) 

at the end. The quark contribution cannot be calculated just from perturbative QCD because 
of low energy strong interaction effects. Rather we determined the contribution, Aq;||^^j(M|), 
by evaluating the dispersion relation (j46p . The results were shown in Table [7| We found 

Aa£d(^l) = 0.02755 ±0.00019e.p±0.00013,ad,vp±0.000019rad,FSR (202) 
= 0.02755 ± 0.00019exp ± 0.00013rad (203) 
= 0.02755 ± 0.00023, (204) 

if we use the inclusive measurements of R{s) in the interval 1.43 < ^/s < 2 GeV. The corre- 
sponding value of the QED coupling is given by 

a{Ml)-^ = 128.954 ± 0.031. (205) 

If, on the other hand, we were to use the sum of the exclusive data for the various e~^e~ — > hadron 
channels, then the result would become 0.02769 ± 0.00018exp ± 0.00013rad and a(M|)-i = 
128.935 ± 0.030. Table [7| shows the contributions to Aa[|^j(M|) from the different energy 
intervals of the dispersion integral, ()46|) . together with the sum. An alternative view may be 
obtained from the (lower) pie diagrams of Fig. ^| They display the fractions of the total 
contribution and error coming from various energy intervals in the dispersion integral. As 
anticipated, both Table[7|and the pie diagrams, show that the hadronic contributions to a(M|) 
are more weighted to higher s values in the dispersion integral for Aa[^^^^(M|), than those in 
the integral for a^^'^ needed to predict (7 — 2 of the muon. 

The above values of Aa^^^{M'^), and the corresponding values of a^^ at s = M|, are 
compared with other determinations in Fig. 123 The BES data [HH] became available for the 
analyses from |133j onwards. In Table fTUl we compare contributions to the dispersion relation Q 
for Aa;[^j(M|) obtained in this work with those found by Burkhardt and Pietrzyk |100j . Since 
new e~^e~ data became available for the former analysis, the comparison is only meaningful in 
the higher energy intervals. Nevertheless, although the agreement in the size of the contributions 
is good, we see that the latter analysis has considerably larger uncertainties in some energy 
intervals, which explains, in part, the difference in the size of the overall error shown in Fig. 123 

8.2 Imphcations for the global fit to electroweak data 

The value of the QED coupling on the Z pole is an important ingredient in the global flt of all 
the precise electroweak data. The continuous curve in Fig. QUI shows the proflle as a function 
of InrriH obtained in the global analysis if our value of Aahad is used, (whereas the dashed 
shows the proflle that would result from the BPOl |100j determination of the QED coupling). 
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Figure 25: Recent determinations [Ml 11001 1103 llH3 1103 [Ml MB [HOI [IHl [122] of 

Aa|j^^^(M|) (lower scale) with the corresponding value of a(M|)~^ at the Z boson mass shown 
on the upper scale. The last two entries, HMNT(03), are the values obtained in this work, and 
include the recent CMD-2 (re-analysed) data jTO] in the evaluation. 



59 



6 



2- 







1 \ \ ■■■ 

' '■ \ — 


0.02761±0. 
0.02755±0. 


00036 / 
00023 / 

i / 

' / 








Excluded 


— ' — ■ r">'i-'^ 


./ Preliminary 




10 


2 



mn [GeV] 



Figure 26: The profile versus the mass of a Standard Model Higgs boson obtained in a global 
analysis of electroweak data. The solid curve is obtained using the value we found in this work, 
and the dotted curve is obtained using the value in |100j . We thank Martin Griinewald for 
making this plot. 

The measured value of mj has been included in the analysis. When our new determination is 
taken, the fit predicts that a Standard Model Higgs has a mass 

mu = 102 +3^gS GeV (206) 

or rriH < 221 GeV at the 95% confidence level. 



9 Conclusions 

The anomalous magnetic moment of the muon, {g — 2)/2, and the QED coupling at the Z 
boson mass, a(M|), are two important quantities in particle physics. At present, the accuracy 
of the theoretical predictions is limited by the uncertainty of the hadronic vacuum polarization 
contributions. Here we use all the available data on e^e^ hadrons to achieve the best 
presently possible data-driven determination of these contributions. In this way, we obtain a 
Standard Model prediction of the muon anomalous magnetic moment of 

af^ = 0.00116591763(74), (207) 

to be compared with the present experimental value of 

a^^P = 0.0011659203(8), (208) 

which shows a 2Aa difference. As this comparison of the measurement and prediction becomes 
more and more precise, we will obtain an increasingly powerful constraint on physics beyond 
the Standard Model. 



60 



We have also used our optimal compilation of the available e'^e 



hadrons data to predict 



a{Mlr 



1 



128.954 ±0.031. 



(209) 



The accuracy is now 24 x 10~^. This again is an important quantity. It is the most poorly- 
determined of the three parameters which specify the electroweak model. Although significantly 
improved from the error of Burkhardt and Pietrzyk's preliminary result [T72 ] , it is still the 
least accurately determined of the three fundamental parameters of the electroweak theory; 
AGJG^ = 1 X 10-5 and AMz/Mz = 2 x 10"^ 

9.1 Future prospects for reducing the error on g — 2 

We have stressed that the comparison of the measurement and the Standard Model prediction 
of the muon anomalous magnetic moment, = ((7^ — 2)/2, is very important. It provides a 
valuable constraint on, or an indicator of, new physics beyond the Standard Model. From the 
above discussion, we see that the present uncertainties on the measurement and the prediction 
are 8 and 7 x 10"^° respectively. How realistic is it to improve the accuracy in the future? On 
the experimental side, the accuracy is dominated at present by the BNL measurement. We 
can expect a further improvement in the BNL measurement of {g — 2), since the collaboration 
are at present analysing 3.7 billion fi~ events which should give a total relative error of about 
0.8 ppm. As a consequence, the ±8 x 10^^° uncertainty in (j208|) should be improved^'^ to about 
±6 X 10^^°. If the error on the theory prediction can be improved beyond this value then the 
case for another dedicated experiment with even more precision is considerably enhanced. 

The error attributed to the theoretical prediction of is dominated by the uncertainties in 
the computation of the hadronic contribution, a^^'^; in particular on the calculation of aJ^^'d.LO 
and a^^'^'^'^^, which at present have uncertainties of about 6 and 4 xlO"^'' respectively. The 
latter error, on the light-by-light contribution, is generally believed to be able to be improved 
to 2 xlO~^° (25% error); and, optimistically, it is perhaps not hopeless to envisage an even- 
tual accuracy of about 1 xlO^^'^ (10% error), but this would require a breakthrough in the 
understanding of this contribution. We are left to consider how much the error on a^^'^'^'^ could 
be improved. Already we are claiming a 1% accuracy. To reduce the error from the present 
6 xlO^^*^ to 1 xlO^^° is not realistic. However we should note (see, for example, Ref. |173j ) 
there will be progress from all experiments that are measuring R. Indeed, with the improve- 
ments, already in progress or planned, of the BES, CMD-3 + SND at VEPP-2000, BaBar, 
Belle, CLEO-C and KLOE experiments, we may anticipate an eventual accuracy of 0.5% in the 
crucial p domain and 1-2% in the region above 1 GeV. It will be challenging, but not impossible. 
This statement also applies to improving the accuracy of the radiative corrections. 

In this connection, note that the measurements of the radiative return experiments are 
just becoming available. From these experiments we may anticipate low energy data for a 
variety of e^e~ channels, produced via initial state radiation, at the 0-factory DA$NE ^ |U] 
and at the S-factories, BaBar and Belle, see, for example, |174j . For instance, by detecting 
the 7r+7r^7 channel, it may be possible to measure the vital e^e~ tt+tt^ cross section in 

^^See the note added in proof. 
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the threshold region. For the radiative return experiments there is no problem with statistics, 
and the accuracy is at present due to systematics, which come mainly from theory. These new 
experiments are motivating much theoretical work to improve their accuracy. Already, today, 
it is claimed to be 2% in the p region. 

In summary, we may hope for an improvement in accuracy down to about 3 x 10~^° in the 
theoretical prediction of in the foreseeable future, which in turn emphasizes the need for an 
experimental measurement with improved precision. 

Note added in proof 

The BNL muon g — 2 collaboration have just published |178j the results of their analysis 
of the fi~ data which updates their experimental determination of a^. As a result they now 
obtain a new world average 

a^^P = 0.0011659208(6) . (210) 

Comparing this value with our SM prediction of Eq. (|207p we find a 3.3a discrepancy, as 
shown by the HMNT (03) (incl.) error bar in Fig. |2ZI That is the discrepancy is Sa^ = 
(31.7 ±9.5) X 10^1°. 

Also, very recently, the hadronic light-by-light contribution has been recalculated, paying 
particular attention to the matching between the short- and long-distance behaviour |179j . The 
contribution is found to be a^^'^'^'^'^ = (13.6 ±2.5) x 10~^°. In addition Kinoshita and Nio have 
updated the calculation of the a"^ QED contribution and find jl80j 

aj^™ = 116584719.35(1.43) x 10~^\ (211) 

which should be compared with the value ()195|) we have used. If we use these new a^^'^'^~^~^ 
and a|^™ values, then our prediction is given by the HMNT (03b) (incl.) error bar in Fig. |27| 
and corresponds to 

af^ = (11659183.53 ± 6.73) x 10~^° (212) 

in the place of Eq. ()207p . If this prediction is compared with the new BNL result above, then 
there is a discrepancy of 2.7a, that is 6a^ = (24.5 ± 9.0) x 10^^°. 
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Figure 27: The new (world average) experimental value of the muon ((^^j — 2)/2 = given 
in Ref. |178j . compared with the SM prediction as given in the text, HMNT (03), and with 
the value, HMNT (03b), which is obtained using the hadronic light-by- light contribution as 
recently calculated in Ref. |179j and the updated QED contribution given in Ref. |180j . 

Appendix A: Threshold behaviour of 77*^7 and 777 produc- 
tion 

We take the Wess-Zumino-Witten (WZW) local interaction as 

i^wzw = --^cpPF^^F^-", (213) 

where /^r ~ 93 MeV, and P denotes the electrically neutral members, 7r° or rjs, of the SU{3) 
pseudoscalar octet. The cp coefficients are 0^,0 = 1 and c^g = 1/a/3. We may extend the 
multiplet to include the SU{3) singlet, r/i, for which the coefficient is c,,^ = 2\/2/\/3. As usual, 
F^jy is the QED field strength tensor, and is its dual, 

F^u = ^pvpcfF^ , (214) 

where e^j^po- is a totally antisymmetric tensor with eoi23 = 1- 
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A.l: TT^ 27 decay and e~^e 7r^7 

The WZW interaction, ()213|) . is responsible for the 7r° — >• 27 decay. The lowest-order amplitude 
M is 

M = ^ (p^)^ ^215) 



which results in the partial decay width 



2 'H 



- 27) = (216) 

when summed over the polarization of the final state photons. If we take fj^ = (130 ± 5)/^/2 
MeV and m^o = 134.9766 ± 0.0006 MeV IH^, then this gives 



r(7r° ^ 27) = 7.81 ± 0.60 eV, (217) 
which is in good agreement with the experimental value |104j . 

r(7r°^27) =7.7±0.6eV. (218) 

exp. 

The cross section of e~^e~ 7r°7 can be written in terms of the vr^ 27 width as 



a(e~^e 



8a7rr(7r° ^ 27) / m 



2 



3 



. A) = cr^e^e- - A) = ^^^^3 ' ' ) ' ^^19) 

We can further improve the behaviour of the cross section by assuming vector meson dominance: 

avMD(e+e- ^ A) = ^pt(e+e- ^ A) (^^) ■ (220) 

\mi - sj 

We use the equation above in calculating the n^'j contribution from the threshold region in 
Section 13.11 

A. 2: ?7 ^ 27 decay and e+e~ rjj 

If we neglect rjQ-rji mixing and identify r/g as r], then the 77 — > 27 decay is dictated by the WZW 
interaction, 

a ~ 

Cwzw = - ^8 F^,F^\ (221) 

which contains an extra factor of 1/ as compared with the 71^77 coupling term. The calcu- 
lation of the decay rate is exactly analogous to that of tt" decay. The result is 

2 s 
(y. Tfi 

r(r/ ^ 27) = —^—^ (LO ChPT without r/i-r/g mixing). (222) 
1927r^/^ 

Taking = (130 ± 5)/V2 MeV and = 547.30 ± 0.12 MeV [TUl, we obtain 

T{r] 27) = 0.174 ± 0.013 keV (LO ChPT without r^i-r^g mixing), (223) 
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(225) 



which differs from the observed value |lU4j by about a factor of 3, 

T{ri^ 27) l^^p = 0.46 ± 0.04 keV. (224) 

The agreement becomes better when we allow for the mixing between the 77 and rj' states. 
Following Ref. |175j . we define the mixing angle 6p by 

/ ?7 \ / cos 6'p - sin 6'p \ / ?78 \ 
\ Tj' / \ sin^^p cos6p / \ Vi J 
The Lagrangian now becomes 

C = -—— (c^, cos 9p - c^, sin Op) ^F^.F^" - -— (c^, sin 9p + Cr,, cos 9p) v'F^.F^". (226) 

If we take 6p ^ —20° [1.75j|. then the coefficient of the rjFF term is 

c^g cos 6'p — c^^ sin ^P = 1.91 x c^^ = 1.10, (227) 
and the predicted decay width is 

(y 771 

T{7]^2-f) = g4^3 j2 (^^8 COS - Cm sin Op)^ 

~ 0.63 keV (LO ChPT with r/i-r/g mixing). (228) 

We find the residual discrepancy with the observed rate is removed when we introduce the 
higher-order effect, fi 7^ fs fn- In this case, 

^ = i^^osOp -^smOp] 7^F,,F'^^ - ^ (^sinOp + ^ cos Op] 7^'F,,F^\ (229) 

If we take /g ^ 1.3/^, /i ^ 1.1/^, as given by Eqs. (162) and (163) of Ref. [Hni, and 9p ^ -20°, 
then the Lagrangian becomes 

01 ~ 

C ~ -1.60 X c^^^F^^F^\ (230) 
and the predicted decay rate is 



r(r/-.27) = ^^-^ l^^c^gcos^p- ^c^,sin^j^ 

~ 0.45 keV (NLO ChPT with r^i-r/g mixing), (231) 

which is now in excellent agreement with the observed value, ()224|) . 

Similarly to the e~^e^ 7r°7 case, we can use the VMD approach to predict the cross section 
of e'^e'^ —>■ Tj'-f. 

( 771^ \ ^ 

crvMD(e+e" ?77) = o-pt(e+e" 777) „ , (232) 

* m,-^, — s 



where 



apt (e+e- ^ r/7) ^ ^ (^^ cos ^^P - ^ sin Op^ ' (l - ^) ' • (233) 

We take the parametrization ()232|) in calculating the e+e^ —>■ Tj'y cross section near the threshold 
region in Section 13.41 
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Appendix B: Constraints on V ^ decay branching frac- 
tions 




7 



V 

\ \ 
q — p \ X (j 

Figure 28: The V ^ 'ja decay in the VMD approach. 

Here we calculate the V —>■ 70" decay of a vector meson using the Vector Meson Dominance 
(VMD) model. To calculate the amplitude, we have used the VMD Lagrangian |177j 

>CvMD = -If^^F^'" - ^V^uV^" + ^ml^V^V^ - gv^nV^J^. - eJ^A^ - (234) 
4 4 z ^9v 

where J ^ is the electromagnetic current. V^^, is defined by 

V^, = d,V, - d,V^. (235) 

Here describes the neutral vector meson (V" = p,uj,(j),- ■ ■). We take e to be positive. The 
diagram which contributes to the decay is shown in Fig. |2H1 The amplitude Ai is given by 

= -Aig^^^ — {{p ■ q)ga(3 - Pag/3)ef (p)ey(g), (236) 
9v 

where ey and are the polarization vectors of V and the photon, respectively. We have 
assumed that the interaction between the a meson and photon is given by 

C = g^^^aF^^F^'", (237) 

where ga^^^ is a coupling constant. From the amplitude of (|236p we can readily calculate the 
required partial decay width 



If we use the parameters 



= 1019 MeV, = 4.26 MeV, 

I/tt = 14.4, 5(0 ^ 7a) < 0.002, (239) 



and assume mo- = 600 MeV, then the coupling constant (70-77 is constrained to be 

g^^^ < 5.2 X 10-^(MeV-^). (240) 
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This bound gives constraints on B{uj cr'j) and -8(0(1.68) cr'j). From ()238p . the branching 
ratio B(y (J7) is 



For the u decay, we use the parameters 

= 783 MeV, = 8.44 MeV, gl/n = 23.2, = 600 MeV, (242) 

to obtain the constraint 

B{u a-f) < 7.2 X 10"^ (243) 
Similarly, for the 0(1.68) — > (17 decay, we have 

5(0(1.68) ^ (T7) < 3.5 X 10"^ (244) 

using the parameters 

m^(i.68) = 1680 MeV, F^(i.68) = 150 MeV, 4(i.68)/vr = 249, = 600 MeV. (245) 

These constraints are used in Section 
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Unannei 


bxperiments with Keierences 




OLYA [IHlinill8^, OLYA-TOF JJj, NA7 [20^ OLYA and CMD [21, 
122, DM1 [23, DM2 [211, BCF [23123, MEA [23 I2H1, ORSAY- 
ACO [23, CMD-2 [inillllllll 


7r°7 
VI 


SND [HilISS 

SND [S21IS3, CMD-2 iMISSi 36j 




ND [22], DM1 [SH, DM2 [HH], CMD-2 [El [TBlinili^ , SND gOlE^, 

CMD [i^ 


K+K^ 


MEA [2Ij, OLYA [13, BCF [23, DM1 [Sj, DM2 [13113, CMD (22j, 
CMD-2 P', SND TTl 




DM1 48J, CMD-2 ^[111113, SND [17| 




MSN DM2 OLYA [S2], CMD-2 [23, SND [M|, ORSAY- 
ACO [S3, 772 [S3, MEA [SZl 




ND and ARGUS 22, DM2 CMD-2 [S3 IS3, SND [Si 103, 

ND [nn 




ND 123, M3N [S3, CMD [03, DM1 [03103, DM2 [Slj, OLYA [03, 
772 [03, CMD-2 p, 67, "03, SND [SI], ORSAY-ACO [S3 


7r^7r~7r"'"7r~7r° 


MEA [Ff], MSN [SDj, CMD [221103 , 772 [56j 




MSN 




DM2 [S3, CMD-2 [69J, DM1 70J 


TT ^ TT ^ TT ^ TT ~ TT TT ~ 

7r+7r^7r'^7r°7r'^7r*^ 


MSN [S3, CMD DM1 DM2 [13 

MSN [S3, CMD [03, DM2 [13, 772 [S3, MEA [SI] 

isospin-related 




DM2 [13, CMD-2 [69J 




DM2 74, 75 


KinK 


DM1 [Z3, DM2 [Z11IZ3 


KIX 


DM1 [77| 




DM2 [73 


pp 


FENICE [TH'.'Tg, DM2 [HOI EH, DM1 [H3 


nn 


FENICE [78, ,83 


incl. (< 2 GeV) 


772 |84i, MEA [H3, MSN "H^, BARYON-ANTIBARYON IHTl 


incl. (> 2 GeV) 


EES [EH IH3, Crystal Ball 1201 HH [13, LENA j93j, MD-1 [^3, 
DASP [23, CLEO [23, CUSB [23, DHHM pi 



Table 1: Experiments and references for the e~^e~ data sets for the different exclusive and the 
inclusive channels as used in this analysis. The recent re-analysis from CMD-2 lOJ supersedes 
their previously published data for 7r^7r~ [TT], 7r'''7r~7r° [13 and K^K^ [Tlj . 
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Experiment 


Procs. 


Norm. 


|cos6'cut| 


Type 


Cvp(a 


spacel.) 


NA7 20_^ 






— 


B 


1.000 


OLYA [iniiniiiHiETiEaiiHi 


TT^TT^, KK 


ee + /i/i 


< 0.71 


D 


0.998 


- 0.993 


[SI ESI 


4:71 






D 


0.995 


- 0.993 


CMD laiEs 


TT~^TT~ , KK 


ee + /i/i 


< 0.60 


D 


0.999 


- 0.994 




3tt, 4tt 






D 


0.996 


- 0.994 


OLYA-TOF [ini 




ee + /i/i 


< 0.24 


D 


0.999 


- 0.998 


MEA ^ 


TT + TT^, KK 


ee 


< 0.77 


D 


0.992 




[2H! 




jjLjJ, 


— 


B 


1.000 




EH 


4:71 


ee 


< 0.77 


D 


0.993 


- 0.992 


DMi [23 El EH] 




ee 


< 0.50 


D 


0.998 


- 0.994 


[^inainii 


Stt, 47r 






D 


0.998 


- 0.994 


DM2 [23113110! 




jJLjJ, 


— 


B 


1.000 




[SHUSH 


Svr, 47r 


ee 


unknown 


— 


no corr. appl. 


SND [SH E21 EI 


7r°7, JsTfs: 


ee 


(< 0.89) 


A 


0.974 


-0.967 


[innmisii 


Svr, 47r 






A 


0.973 


- 0.963 


CMD-2 [HIISII 


KK 


ee 


(< 0.64) 


A 


0.968 


-0.967 


^^^^^1 ^^^rn ^^^^^1 ^^^^^1 ^^^T^i ^^^^^1 

[131 153| l39| 1531 Ib7| IbSj 


Svr, 47r 






A 


0.972 


- 0.963 


772 [HH 


i? 


ee 


< 0.64 


E 


0.992 


-0.991 


DASP [95 


i? 


ee 


< 0.71 


E 


0.985 


DHHM [98J 


R 


ee 


< 0.70 


D 


0.990 


- 0.989 


EES [88, 89j 


R 


ee 


(< 0.55) 


B 


1.000 


Crystal Ball [SOI IHH 112] 


R 


ee 




B 


1.000 


LENA [9"3^ 


R 


ee 




B 


1.000 


CLEO 1^ 


R 


ee 


various 


B 


1.000 



Table 2: Information about vacuum polarization correction factors for different data sets as 
explained in the text. The letters A, B, D, E indicate that the correction factor is given by 
ffTUI) . ([TT|) . (fTS|) . (fTB|) respectively. The tt+tt" and most recent tt+tt^tt^ and K^K^ data from 
CMD-2, as well as the R measurements from BES are given as undressed quantities and are 
already corrected for vacuum polarization effects. According to their publications also the R 
data from CLEO, LENA and Crystal Ball have leptonic and hadronic VP corrections applied 
both in the Bhabha and the hadronic cross sections. The correction factors of type A, D and E 
displayed in the last column are obtained using a{—s) as an approximation to a{s). However 
in the actual analysis we evaluate the corrections using a{s), see Fig. [H 
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channel 




TT + TT Tl^ 


TT^TT TT + TT 


TT+TT vr^vr^ 




AaJ X 10^° 


+ 1.77 


-0.68 


-0.10 


-0.28 




A(A«had(M|))-P X 10^ 


+0.06 


-0.07 


-0.02 


-0.05 




channel 


K+K~ 


mi 




inch (< 2 GeV) 


inch (> 2 GeV) 


Aa;;p X 10^° 


-1.05 


-0.17 


-0.16 


-0.54 


-0.07 


A(Aahad(M|))-P X 10^ 


-0.14 


-0.02 


-0.01 


-0.18 


-0.54 



Table 3: Shifts of the contributions to and Aahad(Af|) from the different channels due to the 
application of the appropriate vacuum polarization corrections to the various data sets. The 
values Aa^P are derived as the difference of calculated with and without VP corrections. 



channel 


data range 


5 (MeV) 


Xmin/d-O.f. 


range used 




Aa^ 


w/o fit 




0.32-3 


3.5 


1.07 


0.32 - 1.425 


502.76 


5.01 


500.10 




0.483 - 2.4 


20,0.6,0.6 
20,0.2,0.2 


2.11 
1.44 


0.66 - 1.425 
0.66 - 1.425 


46.05 
46.42 


0.63 
0.76 


46.54 
47.38 




0.765 - 2.245 


11 


2.00 


0.765 - 1.432 


6.18 


0.23 


5.70 




0.915 - 2.4 


10 


1.28 


0.915-1.438 


9.89 


0.57 


9.44 


K+K- 


1.009 - 2.1 


5,0.6 


1.00 


1.009 - 1.421 


21.58 


0.76 


21.31 


mi 


1.004- 2.14 


10,0.1 


0.86 


1.004 - 1.442 


13.16 


0.16 


13.11 


inclusive 


1.432 - 3.035 
2 - 11.09 


20 
20 


0.28 
0.74 


1.432 - 2.05 
2 - 11.09 


32.95 
42.02 


2.58 
1.14 


31.99 
41.51 



Table 4: Details of the clustering and fit for the dominant channels as described in the text. 
The values of and its error have been multiplied by 10^° and energy ranges are given in 
GeV. For the tt+tt^tt^ channel the bands of clustered data for uo and displayed in Fig. IHl 
were obtained using a clustering size of 0.6 MeV, which leads to a slightly worse Xmin' but a 
better eyeball fit, than for the 0.2 MeV clustering. For the numerics we have used the 0.2 MeV 
clustering size. The differences are small. 
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channel 


1 • / -1 

inclusive (J 

had T,0 


.43,2 GeV) 

Afthad (M|) 


1 • / -1 

exclusive (. 

had T,0 


-.43,2 GeV) 
Aahad(M|) 


7r"7 (ChPT) 


0.13 ± 0.01 


0.00 ± 0.00 


0.13 ± 0.01 


0.00 ± 0.00 


TT 7 (data) 


4.50 ± 0.15 


0.36 ± 0.01 


4.50 ± 0.15 


0.36 ± 0.01 


TT+TT (ChPT) 


2.36 ± 0.05 


0.04 ± 0.00 


2.36 ± 0.05 


0.04 ± 0.00 


TT^TT (data) 


502.78 ± 5.02 


34.39 ± 0.29 


503.38 ± 5.02 


34.59 ± 0.29 


TT+TT TT^ (ChPT) 


0.01 ± 0.00 


0.00 ± 0.00 


0.01 ± 0.00 


0.00 ± 0.00 


-1- — / 1 J \ 

TT^TT TT (data) 


46.43 ± 0.90 


4.33 ± 0.08 


47.04 ± 0.90 


4.52 ± 0.08 


777 (ChPT) 


0.00 ± 0.00 


0.00 ± 0.00 


0.00 ± 0.00 


0.00 ± 0.00 


777 (data) 


0.73 ± 0.03 


0.09 ± 0.00 


0.73 ± 0.03 


0.09 ± 0.00 




21.62 ± 0.76 


3.01 ± 0.11 


22.35 ± 0.77 


3.23 ± 0.11 




13.16 ± 0.31 


1.76 ± 0.04 


13.30 ± 0.32 


1.80 ± 0.04 


27r+27r" 


6.16 ± 0.32 


1.27 ± 0.07 


14.77 ± 0.76 


4.04 ± 0.21 


-i- — r\ n 


9.71 ± 0.63 


1.86 ± 0.12 


20.55 ± 1.22 


5.51 ± 0.35 


^ _i_<r\ n 

27r+27r tt'^ 


0.26 ± 0.04 


0.06 ± 0.01 


2.85 ± 0.25 


0.99 ± 0.09 


^ — D 


0.09 ± 0.09 


0.02 ± 0.02 


1.19 ± 0.33 


0.41 ± 0.10 




0.00 ± 0.00 


0.00 ± 0.00 


0.22 ± 0.02 


0.09 ± 0.01 


27r+27r"27r° 


0.12 ± 0.03 


0.03 ± 0.01 


3.32 ± 0.29 


1.22 ± 0.11 


TT^TT 47r" (isospm) 


0.00 ± 0.00 


0.00 ± 0.00 


0.12 ± 0.12 


0.05 ± 0.05 


K^K tt" 


0.00 ± 0.00 


0.00 ± 0.00 


0.29 ± 0.07 


0.10 ± 0.03 


K^i^^^TT^ (isospm) 


0.00 ± 0.00 


0.00 ± 0.00 


0.29 ± 0.07 


0.10 ± 0.03 




0.05 ± 0.02 


0.01 ± 0.00 


-I r\ /~\ 1 /~\ -I 

1.00 ± 0.11 


0.33 ± 0.04 


ziz T^-\- /• • \ 

K^TT+K^ (isospm) 


0.05 ± 0.02 


0.01 ± 0.00 


-I r\ /~\ 1 -1 -1 

1.00 ± 0.11 


0.33 ± 0.04 


KKttit (isospm) 


0.00 ± 0.00 


0.00 ± 0.00 


3.63 ± 1.34 


1.33 ± 0.48 


UJ[—>- TT 7)7r 


0.64 ± 0.02 


0.12 ± 0.00 


0.83 ± 0.03 


0.17 ± 0.01 


/ n \ 4- 

a;(^ TT 7)7r^7r 


0.01 ± 0.00 


0.00 ± 0.00 


0.07 ± 0.01 


0.02 ± 0.00 


/ \ 4- 

77(— *• TT 7)7r^7r 


0.07 ± 0.01 


0.02 ± 0.00 


0.49 ± 0.07 


0.15 ± 0.02 


(/)(^ unaccounted) 


0.06 ± 0.06 


0.01 ± 0.01 


0.06 ± 0.06 


0.01 ± 0.01 


pp 


n nn -i- n nn 
U.UU It u.uu 


n nn -i- n nn 
u.uu It u.uu 


n n/1 -I- n ni 

U.U4 It U.Ui 


n no -1- n nn 
u.u/ It u.uu 


nn 


0.00 ±0.00 


0.00 ±0.00 


0.07 ±0.02 


0.03 ±0.01 




7.30 ± 0.43 


8.90 ±0.51 


7.30 ± 0.43 


8.90 ±0.51 


T(15-65) 


0.10 ±0.00 


1.16 ±0.04 


0.10 ±0.00 


1.16 ±0.04 


inclusive R 


73.96 ± 2.68 


92.75 ± 1.74 


42.05 ± 1.14 


81.97 ± 1.53 


pQCD 


2.11 ±0.00 


125.32 ±0.15 


2.11 ±0.00 


125.32 ±0.15 


sum 


692.38 ±5.88 


275.52 ±1.85 


696.15 ±5.68 


276.90 ± 1.77 



Table 5: Contributions to the dispersion relations and (0) from the individual channels. 
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A/i (GeV) 


comment 


^had,LO y XQ^Q 


Aahad(M|) X 10^ 


0.32-1.43 




502.78 ± 5.02 


34.39 ± 0.29 


(0.32-1.43 


'old' CMD-2 


492.66 ± 4.93 


33.65 ± 0.28) 


0.32-2 




503.38 ± 5.02 


34.59 ± 0.29 


0-0.32 


ChPT 


2.36 ± 0.05 


0.04 ± 0.00 



Table 6: tt+tt contributions to aJ^a^d-LO AQ;had(Af|) from various energy intervals. The 
entries in brackets give the contributions obtained using the CMD-2 data before re-analysis. 



energy range (GeV) 


comments 


^had.LO X iQio 


Atthad (M|) X 10^ 


m,r^0.32 


ChPT 


2.36 ±0.05 


0.04 ±0.00 


0.32-1.43 


excl. only 


606.55 ±5.22 


47.34 ±0.35 


1.43-2 


inch only 


31.91 ±2.42 


10.78 ±0.81 




(excl. only 


35.68 ± 1.71 


12.17 ±0.59) 


2-11.09 


inch only 


42.05 ± 1.14 


81.97 ± 1.53 


J/ip and ip' 


narrow width 


7.30 ± 0.43 


8.90 ±0.51 


T{1S-QS) 


narrow width 


0.10 ±0.00 


1.16 ±0.04 


11.09-00 


pQCD 


2.11 ±0.00 


125.32 ±0.15 


Sum of all 


inch 1.43-2 
(excl. 1.43-2 


692.38 ±5.88 
696.15 ±5.68 


275.52 ± 1.85 
276.90 ± 1.77) 



Table 7: A breakdown of the contributions to different intervals of the dispersion integrals for 
jjhad,LO ^j^^ Aahad(^l)- The alternative numbers for the interval 1.43 < ^/s < 2 GeV corre- 
spond to using data for either the sum of the exclusive channels or the inclusive measurements, 
see Fig. IH 
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(a) Breakdown of contributions to l.h.s. of sum rules 



energy range (GeV) 


contribution (m = 2, n = 0) 


contribution {m — n — 0) 


2m^ - 0.32 (ChPT) 


0.00 ± 0.00 


0.00 ± 0.00 


0.32 - 1.43 (excl) 


3.92 ± 0.03 


4.49 ± 0.04 


1.43 - 2.00 (excl) 


3.02 ± 0.26 


4.93 ± 0.43 


1.43 - 2.00 (incl) 


2.48 ± 0.19 


4.03 ± 0.30 


2.00 - 3.73 (incl) 


3.94 ± 0.14 


22.56 ± 0.70 


sum (excl) 


10.87 ± 0.30 


31.98 ± 0.82 


sum (incl) 


10.34 ± 0.24 


31.08 ± 0.76 



(b) Breakdown of contributions to r.h.s. of sum rules 



origin 


contribution (m = 2, n = 0) 


contribution (m = n = 0) 


massless QCD 
correction from finite 
quark and gluon condensates 


10.31 ± 0.05 
-0.03 ± 0.02 
0.03 ± 0.02 


30.43 ± 0.11 
-0.03 ±0.02 
0.00 ± 0.00 


prediction from QCD (total) 


10.30 ± 0.06 


30.40 ± 0.12 



Table 8: The breakdown of the sum rules for ^/so = 3.7 GeV, for the choices m = 2,n = and 
m — n — 0. The contributions to the left-hand-side (data) are shown in the upper table, and 
the QCD contributions are given in the lower table. 
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channel 


this work (^s < l.SGeV) 


DEHZ 03 (yi < 1.8GeV) 


difference 


■n+TT- (ChPT) 


2.36 ±0.05 (< 0.32GeV) 


58.04 (± 2.06) 


(<0.5GcV) 




TT+TT^ (data) 


503.24 ±5.02(> 0.32GeV) 


450.16 (± 5.14) 


(>0.5GeV) 




TT+TT^ (total) 


505.60 ± 5.02 


508.20 ± 5.53 




-2.60 


T! 


0.13 ± 0.01 (ChPT) 


0.93 








4.50 ± 0.15 (data) 


+37.96 X 0.0889 


{u 7r07) 








+35.71 X 0.00124 


(0 ^ ^07) 




777 


0.01 ± 0.00 (ChPT) 


+37.96 X 0.0007 


(cj 






0.73 ± 0.03 (data) 


+35.71 X 0.01299 


{<P Vl) 




7r'^7 + ?77 


5.36 ± 0.15 


= 4.84 ± 0.18 




+0.52 


— T n 


0.01 (± 0.00) (ChPT) 


37.96 X 0.9104 {uj - 


-> Tr+TT^Tr*^) 








+4.20 (0.81 < 


Vi< 1.00) 






+46.97 (± 0.90) (data) 


+35.71 x 0.155 (0- 










+2.45 (1.055 < ^ 


/s < 1.800) 






= 46.98 ± 0.90 


= 46.74 ± 1.09 




+0.24 




22.29 ± 0.76 


4.63 + 35.71 X 0.492(0^ 


K+K-) 








= 22.20 ± 0.59 




+0.09 




1 on 1 f\ 00 


0.94 + 35.71 X 0.337(0- 










12.97 ± 0.31 




+0.32 


If /, ()_, 

(P(7^ Sn,2K,Tr j,ri-f) 


U.Ud ± (J. Do 


35.71 X 0.002(0 74 37r, 2/4', 7r"7, 777) 








= 0.07 ± 0.00 




—0.01 


+ — 00 


18.34 ± 1.08 


16.76 ± 1.33 


+1.58 


/ \ 

cj(-^ 71*^7)71*^ 


0.82 ± 0.03 


0.63 ± 0.10 


+0.19 


4- — + — 


13.63 ± 0.70 


14.21 ± 0.90 


-0.58 


+ — + — 
TT^TT TT^TT TT 


2.05 ± 0.18 


2.09 ± 0.43 


-0.04 


4- — 


0.85 ± 0.30 


1.29 ± 0.22 (isospin, 77) 


-0.44 


7 n — \ — X 

w(— » 7r"7)7r^7r 


0.06 ± 0.01 


0.08 ± 0.01 


-0.02 


TT"*" TT ~ TT"*" TT ~ 7r+ TT ~ 


0.07 ± 0.01 


0.10 ± 0.10 




-0.03 


-4- — 4- — n n 

TT^TT TT^TT TT TT 


1.96 ± 0.18 


1.41 ± 0.30 




+0.55 




0.07 ± 0.07 (isospin, r) 


0.06 ± 0.06 (isospin, t) 




+0.01 


sum from Gtt 


2.11 ± 0.19 


1.57 ± 0.34 




+0.54 




0.43 ± 0.07 


0.54 ± 0.07 


-0.11 




0.85 ± 0.09 








KI-kK 


0.85 ± 0.09 (isospin) 








K^ttK + K'^ttK 

J. i / I J v 1 J v r / 1 J v 


1.71 ± 0.19 


1.84 ± 0.24 




—0.13 


K+K-tt" 


0.18 ± 0.05 










0.18 ± 0.05 (isospin) 










0.36 ± 0.11 


0.60 ± 0.20 




-0.24 


KK'K'K 


2.38 ± 0.98 (isospin) 


2.22 ± 1.02 


+0.16 


total (Vi < l.SGeV) 


636.29 ± 5.43 


636.85 ± 6.08 


-0.56 



Table 9: The contributions of the individual e^e channels, up to y/s = 1.8 GeV, to dispersion 
relation (gH) for a^^'^'^^ (xlO^°) that were obtained in this analysis and in the DEHZ03 study [Sj. 
The last column shows the difference. "Isospin" denotes channels for which no data exist, and 
for which isospin relations or bounds are used. We have divided the DEHZ u contribution into 
the respective channels according to their branching fractions |lU4j , with their sum normalized 
to unity. 
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energy range (GeV) 


HMNT 03 


BP 01 


1.05-2.0 


16.34 ± 0.82 (excl+incl) 
(5.56 ± 0.13 (1.05-1.43 GeV, excl)) 
(10.78 ± 0.81 (1.43-2.0 GeV, incl)) 


15.6 ± 2.3 (excl) 


2.0-5.0 


38.13 ±1.10 (incl) 


38.1 ±2.2 (incl) 


5.0-7.0 


18.52 ± 0.64 (incl) 


18.3 ± 1.1 (incl) 


7.0-12 


30.16 ± 0.61 (incl+pQCD) 
( 25.32 ± 0.61 (7.0-11 GeV, incl)) 
( 4.84 ± 0.02 (11-12 GeV, pQCD)) 


30.4 ± 0.4 (incl) 


12-00 


120.48 ±0.13 (pQCD) 


120.3 ±0.2 (pQCD) 



Table 10: Comparison of the contributions to AQ;had(j^l) x 10^ with the analysis of BP 01 |100j . 



83 



